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SUMMARY 


This report contains a description of the Gordian Refinery Simula- 
tion Model. This is a Fortran computing program for predicting the flow 
streams and material, energy, and economic balances of a typical petro- 
leum refinery, with particular emphasis on production of aviation turbine 
fuel of varying end point and hydrogen content specifications. The program 
has provision for shale and coal oil in addition to petroluem crudes; 26 
crude assays are supplied, and the user may add up to 9 more or update the 
existing assays. The program contains relationships for predicting and 
blending the jet fuel freezing point and smoke point as well as internal 
stream and product hydrogen, nitrogen, and sulfur contents and other pro- 
perties. A case study feature permits dependent cases to be run for 
parametric or optimization studies by input of only those variables which 
are changed from the base case. 

This report is the first of three volumes. It is independent and has 
sufficient detail for the needs of most readers in interpreting the program 
features and output data. Two subsequent volumes contain the mathematics 
and data base details and the program listing. 
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COMPUTER MODEL FOR REFINERY OPERATIONS WITH EMPHASIS 
ON 7ET FUEL PRODUCTION: VOLUME I. PROGRAM DESCRIPTION 


1.0 INTRODUCTION 

Major price increases and the impending shortage of petroleum 
reserves with respect to increasing product demand has brought about a 
serious examination of possible changes in jet fuel composition. Speci- 
fication aviation turbine fuel (ASTM D-T 655) is produced from mid-dis- 
tillate petroleum fractions, which compete with ever growing demands for 
diesel, fuel oil and petrochemical feedstocks. Increased distillate production 
from present crudes is feasible, but conversion of gas oils and residuals 
increases the aromatic content of the mid-distillate pool. Moreover, 
promising alternate crude sources, such as shale oil, tar sands, and 
coal liquids yield distillates also with increased aromatic, nitrogen 
and sulfur contents. Special processing would be required to produce 
present specification aviation turbine fuel from these sources. 

This view of the future has stimulated a reexamination of the 
optimum combination of jet fuel specifications, with respect to the 
refinery processing, the supply distribution system, the aircraft fuel 
system, and the fuel combustion qualities. The goals of current studies 
are assessing the suitability of jet fuels produced from cracked petroleum 
and alternate crude sources and developing a data base which will allow 
optimization of future fuel characteristics. Future aviation turbine 
fuel specifications must represent a trade-off between energy and cost 
efficiency of manufacture and aircraft and engine design and performance. 

This report deals with the refinery portion of the overall program. 

In order to have a systematic way of determining the energy efficiency 
of the production of various product slates Involving different crude 
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sources and different processing schemes, the Lewis Research Center of 
NASA has supported the development of this computer model for petroleum 
refinery operation. The primary objectives of this model are: 

1. the flexibility to configure a refinery involving any or all 
of the process units commonly employed in the production of gasoline, 
jet fuels, and mid-distillates; 

2. the ability to produce jet fuel blends of varying end-point' 
specification and varying specified hydrogen content as part of the 
total slate of products; 

3. the ability to handle synthetic crudes (shale and coal derived) 
with varying severities of hydroprocessing; 

4. the determination of overall refinery energy efficiency; 

5. the determination of sulfur, nitrogen, and hydrogen material 
balances for each process unit and for the overall refinery; and 

6. the capability of carrying out economic calculations. 

The Gordian Refinery Simulation Model, presented herein, has all 
the above capabilities. This report, Volume I of three volumes, contains 
the detailed information necessary to run the program, and a brief 
discussion of the technical basis of the model and its limitations in 
the present state. Volume II (NASA CR-1 35334) contains details of the 
mathematical relations and data bases used in the program, and Volume 
III (NASA CR-1 35335) contains programming documentation. The complete 
documentation and program tape are available through the Computing 
Software and Management Information Office (COSMIC) under the number 
LEW- 13047. 
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2.0 PROGRAM DESCRIPTION AND LOGIC 


2.1 General Description and Applications 

The Gordian Refinery Simulation Model is a Fortran computer program 
which is capable of predicting the production volumes of petroleum 
refinery products, with particular emphasis on aircraft turbine fuel 
blends and their key properties. The program also provides the calculation 
of capital and operating costs for the refinery and its margin of profitability. 
This model is noteworthy in that it includes provision for the processing 
of synthetic crude oils from oil shales and coal liquefi cation processes. 

The program provides .highly detailed blending computations for alternative 
jet fuel blends of varying end-point specification. Volumetric and 
weight balances are output along with complete sulfur, nitrogen and 
hydrogen balances. Energy usage is reported in terms of net fuel power 
and steam consumption along with the hydrogen usage. Results are 
reported in SI metric and English units. 

In order to use the program, it is necessary to specify the type 
and volume of each crude oil which makes up the feed to the refinery as 
well as the capacities of the specific refinery processing units which 
comprise the refinery configuration. The complete slate of refinery 
products is then calculated by the program, with particular emphasis on 
accurately predicting the volumes and properties of the jet fuel, distillate 
and residual fuel oil blends produced. The economic calculations require 
the input of crude oil prices, refinery product realizations, electricity 
cost and investment carrying charge. 
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While the model simulates the operation of a refinery, it does not 
optimize It with respect to maximum dollar profit In the way that a 
linear programming model does. The non-linear relationships which have 
been programmed (e.g. for turbine fuel smoke point and freezing point 
prediction) allow for greater accuracy than may be obtained in calculating 
the product blends normally obtainable with the conventional linear 
programing approach. An additional advantage is that less training is 
required and less data Input and program usage complexities are encountered 
in using a Fortran program as opposed to a complex linear programming 

system. 

Optimization and parametric type studies in general, however, may 
be made by using the case study feature which has been incorporated Into 
the model. This feature is used by establishing a base case and then by 
specifying only the changes to the base case which are applicable to the 
case study sequence, thus minimizing input requirements and permitting 
large case study sequences to be run without undue user or computer 

time. 

Two major data base subroutines are contained within the program, a 
Crude Assay subroutine and a Refinery Process Yield subroutine. The 
Crude Assay subroutine is comprised of the yields and properties of 
distillation products from a wide range of crude oils. The Refinery 
Process Yield subroutine contains typical process unit yields, energy 
requirements and stream quality factors for a variety of refinery processing 
units. In addition, an economic data base is contained in the model 
which provides the fixed and variable components of operating costs and 
investment costs for individual refinery process units. 
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Five potential categories of fuels production are considered by the 
Ref i nery Model : 

1. Kerosene/Aviation turbine fuel 

2. Diesel fuel 

3. Direct gas oil sale 

4. Distillate Fuel 

5. Residual Fuel 

6. Gasoline 

The volumes of kerosene, diesel fuel and gas oil to be produced (if 
any) are specified in the input. Any excess kerosene or diesel is 
routed to turbine or distillate fuel while the remaining gas oil is 
either processed or blended directly to residual fuel. 

An aviation turbine fuel blending subroutine calculates product 
blends at end-points of 525’F and 650«F, and any specified intermediate 
end-point. In addition, the user has the option of specifying the 
required hydrogen content for a given end-point turbine fuel, and the 
model will calculate the maximum quantity of fuel that can be produced 
to meet specification. Other properties of the kerosene/turbine fuel 
pool that are calculated are gravity (“API), sulfur content, nitrogen 
content, viscosity, freeze point, smoke point, paraffins, naphthenes, 
aromatics and heat of combustion. 

The distillate and residual fuel oils usually comprise the bulk of 
the fuel oil production. The model calculates the following for both ; 
the distillate fuel and residual fuel oil pools: 

1. Total pool volumes 

2. Weight percent sulfur 

3. Weight percent nitrogen 

4. Weight percent hydrogen 
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5. API gravity 

6. Viscosity 

If desired, sulfur contents may be specified for distillate and 
residual fuel oil blends. The model will then calculate the maximum 
fuel oil volumes available at the specified sulfur contents, along with 
the other fuel oil properties. In order to make particular sulfur 
specifications, it is usually necessary to reject blending components 
from the total fuel oil pools. The model calculates the volume, sulfur 
content, and the other fuel oil properties of the "rejected" portions of 
the distillate and residual fuel oil pools. 

If desired, the production volume and sulfur content of one or two 
distillate fuel oil blends may be specified. The program will calculate 
the composition and properties of these blends and will blend the remaining 
excess distillate components into the residual fuel oil pool. One or 
two residual fuel oil blends of fixed volume and sulfur content may also 
be specified. The composition and properties of these residual fuel oil 
blends (which now contain excess middle distillate components) will be 
calculated along with the volume and properties of the remaining residual 
fuel oil pool. 

Use of the above blend specification options affords considerable 
flexibility in assessing the production capabilities in a given refining 
situation for fuels ranging from aircraft turbine level to residual fuel 
oil. 


2.2 Technical Basis for Model 

A brief summary of the technical background of the computer model 
is presented in this section. This information is given in more detail 
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in Volume II (NASA CR-1 35334). 

The model has storage capability for up to 35 crude assays. At 
present, the program contains within a subroutine 26 assays, including 
three shale oils and one coal syncrude. The petroleum crude assays were 
obtained from Gordian Associates, Inc. in-house data and from assays 
published in the Oil and Gas Journal. The shale and coal syncrude assay 
data were obtained from references 1 and 2. In those cases where a 
crude oil property needed for the program was not available from the 
assay data, either a correlation was used to derive the property or an 
estimate was used. Such cases are noted in Volume II. 

Process unit yield data for the petroleum-derived fractions were 
obtained from Gordian in-house data and from ref. 3, and for the syncrude 
fractions from references 1 and 2. The process unit yield data base 
includes the energy requirements for each unit in the form of steam, 
refinery gas, and electrical power. 

Several physical and chemical property blending correlations that 
are used in this model are summarized as follows: 

1. Hydrogen content, if not specified as input or through the 
assay, is calculated from API gravity and mean average boiling point by 
the method described in ref. 4 for petroleum fractions. The same method, 
but with adjustments to the predicted carbon-hydrogen ratio, was used 
for the shale oil and coal oil cuts. 

2. Heats of combustion were calculated from API gravity, Watson, 
characterization factor, and sulfur content by a correlating relation 
also obtained from ref. 4. Heat of combustion is blended by weight. 
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3. Smoke points were calculated from Watson characterization 
factor by equations derived from data in ref# 5. Smoke point is blended 

volumetrically by a reciprocal blending relation. 

4. Freezing point data, if not available through assay data, were 
estimated from a correlation presented in ref. 5. Freezing points of 
product streams from the process units were based upon changes in freezing 
point from the feed streams in order to tie the product freezing points 

as much as possible to the original crude oil assays. Freezing point 
blending of the jet fuel and mid-distillate blends was based on the 
blending index method of ref. 6. 

5. Sulfur and nitrogen content, viscosity, and PNA (paraffins, 
naphthenes and aromatics) are based on crude assay values and on inspection 
tests accompanying process unit yield data. Viscosities are blended as 
shown in Exhibit 6 of this volume; the other properties here are blended 

by volume. 

6. The economic calculation for the specific refinery configuration 
under study includes investment costs for the individual process units, 
off-site investment costs, fixed and variable operating costs. Crude or 
other feedstream costs, utility costs, product prices and rate of return 
on investment must be input for the particular case. 

Some cautions regarding the use of the Gordian Refinery Simulation 
Model in its present form must be noted. Process unit data and correlations 
involving petroleum-derived fuels have a broad base of years of experience 
in the petroleum industry, and the calculations involving the use of 
petroleum crudes can be expected to have good reliability. Processing 
data and correlations involving shale or coal -derived fuels, however, 
are extremely limited at the present time, and the results calculated 
using these syncrudes would be expected to be much less reliable. 
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Nevertheless, it was desired to have a refinery simulation model which 
would be capable of handling these syncrudes, with the intent of continually 
upgrading the data and correlation bases as such information becomes 
available. 

2.3 Review of Model Output 

In order to provide an illustration of the capabilities of the 
model this section provides a review of the complete computer printout- 
for a sample case. A copy of an actual computer output is shown in 
Exhibit 1. The computer print-out reports appear in the order in which 
the calculation proceeds within the model. While all of the report 
writing capabilities of the program are shown, there are options to 
selectively suppress various reports in a manner to be specified by the 
user. 

The first set of reports give individual process unit material 
balances with respect to weight flow rate, sulfur, nitrogen, and hydrogen. 
Each material balance process unit report is followed by an energy 
utilization report showing the consumption of fuel, steam, electricity 
and hydrogen and their percentage of the total fuel oil equivalent (FOE) 
input to the refinery. Please note that while the hydrogen unit balance 
appears in error, this results from the fact that v/ater feed and carbon 
dioxide product has been excluded from the weight balance as normal 
refinery practice. 

The second set of reports refer to the jet fuel blends which may be 
produced by the combination of the specified refinery crude oil feedstocks 
and refinery process unit configuration. These reports have been generated 
to provide the full extent of jet fuel blending information which is 
available from the model. These blend reports are: 
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exact matches (except where sulfur and nitrogen is truly lost through 
the combustion of catalytic coke) since the process unit data base 
values are expressed in the program as a percentage distribution of 
sulfur and nitrogen between process unit products. Hydrogen distribution 
on the other hand, is not handled in this manner, but is calculated from 
an API Technical Data Book correlation (ref. 4) which relates hydrogen 
content to specific gravity and the mean average boiling point of the 
stream, with adjustment factors applied for coal and shale oil fractions. 

The overall hydrogen of 100.9 weight percent appears good, but is deceptive. 
While the results for units which process petroleum fraction is generally 
good, rather significant discrepancies can be observed for the shale and 
coal oil processing units. The reason here is a combination of relatively 
weak process unit data available at the present time to relate yields to 
hydrogen consumption and downstream product qualities, and probable 
deficiencies in the correlation used to predict hydrogen content for the 
shale and coal oil species. 

A summary of process unit operations is next shown which compares 
the input process unit capacity with the actual feedrate. Process unit 
severities are also given. Separate reports are printed which summarize 
the overall refinery energy efficiency and usage and the properties of 
the blended gasoline pool. 

The final section of the computer printout report gives an overall 
economic summary of refinery operations. The profit margin calculations 
shown are illustrative only since actual crude oil prices, product 
realizations, and investment carrying charges depend on the specific 
refining situation. However, the process unit Investment costs and 
operating costs are contained in the program as data base items. These 
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(a) Properties of the total 525°F endpoint pool 

(b) Properties and quantity of 525°F endpoint pool at the specified 
hydrogen content. (12.8 weight percent in this run). 

(c) Properties of the total 650°F endpoint pool 

(d) Properties and quantity of 650°F endpoint pool at the specified 
hydrogen content. 

(e) Properties of the total pool at the specified endpoint 
(600°F in this run) 

(f) Properties and quantity of the specified endpoint blend at the 
specified hydrogen content. 

Separate auxiliary reports are shown which give the properties of 
the total middle distillate pool (which include cracked olefinic stocks 
which are normally diverted from jet fuel), and the residual fuel oil 
pool . 

A summary material balance report shows the disposition of mass, 
sulfur, nitrogen and hydrogen between feed streams and final products*. 
Overall recoveries are given at the bottom of this report. Note that 
the volumetric recovery is reasonable for a refinery of the specified 
configuration. The weight recovery is less than 100.0 of weight percent, 
reflecting the fact that weight balances for the individual process 
units have not been forced through the proration of discrepancies. The 
lower sulfur and nitrogen recoveries are due almost entirely to crude 
unit imbalances which reflect inconsistencies in the basic crude assay 
data. In general, basic data of this type has not been adjusted through 
internal program computations because to do so would cause the program - 
user to lose sight of those data base areas which could be improved. 
Sulfur and nitrogen balances for the downstream units, however, provide 



* The middle distillate blend reported here is the balance remaining 
after the specified jet fuel blendihas been made. 
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were obtained by a general review of today's petroleum refining costs 
and are expressed in current 1977 dollars, with no built-in inflation. 

Crude prices, product realizations, electricity cost and investment 
carrying charge must be input since they are specific to a given situation. 

Additional information provided in the front of the report printout 
show the input data cards as punched. The items appearing under the 
headers &CDATA, &UDATA and &UDATA1 are for illustrative purposes only — 
data input under these cards allows for user directed changes to any of 
the data base values contained in the crude oil assay data base subroutine 
(CBASE) and the process unit data base subroutine (UBASE). This capability 
can prove valuable for simulating specific refining situations. 

2.4 Program Logic 

An overall block diagram describing the logic of the Gordian Refinery 
Simulation Model is shown in Exhibit 2. A brief functional description 
comprising the Gordian Refinery Simulation Model follows: 

MAIN PROGRAM - The main program blends crude oil properties, 
calculates the refinery volumetric balance, calculates the overall 
refinery material balance, issues summary reports and sets up the 
subroutine calling sequence. 

BLOCK DATA - Contains initialization data used throughout the 
program and refinery process unit cost data used by the economic subroutine. 

CONAP I - Converts from API to specific gravity and back as required. 

ISTN - Sets stream number identities used by the main program and 
the report subroutines. 

CBASE - Contains the crude assay data base. 

UBASE - Contains the refinery process unit yield and stream quality 
data base. 

MATCLC - Performs all material balance calculations required for 
the weight, sulfur, nitrogen and hydrogen balance around each of the 
refinery process units. 
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INITB - Initializes default input value for each base case run. 


UNTREP - The report writer subroutine for the individual process 
unit material balances. 

UNTENR - The report writer subroutine for energy usage consumed by 
each refinery process unit. 

ECON - Does the refinery economics calculation and outputs summary 
reports . 

HYDCAL - Calculates the hydrogen content of every refinery process 
stream as a function of specific gravity and mean average boiling point — 
unless a value has been pre-specified in UBASE or input. 

INITV - Initializes intermediate program variables. 

BLDARY - Sets up blending stream proprieties for each stream used 
in jet fuel, middle distillate and residual fuel blends. 

BLENDF - Contains the calculation sequence for each of the jet fuel 
blend combinations and the middle distillate and residual fuel oil 
blends. 

SS0RT1- Sorts jet fuel blending components according to their 
hydrogen content ranking. 

BLDREP - The report writer subroutine for turbine fuel, middle 
distillate and residual fuel oil blends. 

BPROP - Calculates blend properties and contains freezing point and 
smoke point blending correlations. 

SSORT - Sorts blend components by sulfur content ranking - used for 
calculating middle distillate and residual fuel oil blend combinations. 

BMAX - This subroutine calculates the maximum quantity of jet fuel 
blend production of a specified hydrogen content, or maximum middle 
distillate or residual fuel of a specified sulfur content. 

A more detailed and comprehensive description of the programming 
aspects of the Gordian Refinery Simulation Model is given in Volume II 
(NASA CR-1 35334) and Volume III (NASA CR-1 35335). 


2.5 Summary of Model Options and Flexibility 

Any combination of crude oils represented in the Crude Assay 
subroutine may be selected by specifying their code number and the 
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charge volume of each crude; including shale oil and coal liquefaction 
products. If the feed to the refinery contains crude oils which are not 
represented in the Crude Assay subroutine, then the yields and properties 
of these crudes may be input directly into the model along with their 
volumes. If Individual refinery charge streams, such as purchased gas 
oil, distillate cutter stocks, and naphthas are a part of the refinery 
charge, these may also be specified by giving their feed volumes and 
stream properties. 

The volumes of up to four separate "imported" fuel oil blending 
components (kerosene, diesel, gas oil, and residuum) may be specified 
along with their properties in order to simulate the refinery operation. 
These components enter directly into fuel oil blending without being 
mixed with other streams. 

The block flow diagram of the refinery processing units is shown as 
Exhibit 3. The following process units may be specified when defining 
the refinery process unit configuration: 

(1) Atmospheric crude distillation (petroleum-based, shale oil, 
coal oil) 

(2) Vacuum crude oil distillation (petroleum-based, shale oil, 
coal oil) 

(3) Catalytic gas oil cracking 

(4) Thermal gas oil cracking 

(5) Middle distillate desulfurizer 

/ 

(6) Kerosene Hydrotreater (separate units for petroleum stocks, 
coal based and shale oil stocks) 

(7) Fluid coking 

(8) Vacuum bottoms visbreaking 

(9) Middle distillate and gas oil hydrocracking (separate gas oil 
hydrocracking units for petroleum stocks, coal -based and shale 
oil stocks) 
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(10) Catalytic naphtha reforming 

(11) Alkylation 

(12) Polymerization 

(13) Butane isomerization 

(14) Hydrogen production 

Variables which serve to describe internal refinery operations are 
also a part of the input. The catalytic cracker severity (percentage 
conversion), the naphtha reformer severity (octane), the type of petroleum- 
based gas oil hydrocracking operation (gasoline versus distillate production, 
and the severity of shale oil and coal oil hydrocracking and hydrotreating 
operation must be specified either directly (or will be by default) if 
these units are included in the overall refinery configuration. In a 
complex refinery where internal streams such as diesel oil and coker and 
vacuum gas oil have alternative dispositions, the fraction of each such 
stream going to alternative refinery process units may be specified. If 
fractional dispositions are not specified, the default is to prorate 
internal refinery stream to alternative down stream units by the ratio 
of their stated capacities. 

In some instances It may be desirable to replace the typical refinery 
process yields contained in the Refinery Process Yield subroutine with 
yield and quality numbers which more exactly represent a particular 
refinery. This could occur if a particular refinery was atypical in 
some respects or if a rather close match of unusual refinery operations 
was desired. For this reason, the option exists to override any of the 
data base values contained in the Refinery Process Yield subroutine by 
specifying them in the input. 
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With regard to the refinery process unit capacities, it would be 
unusual for a specified unit throughput to exactly match the volume of 
an intermediate refinery stream. If the unit is undersized, the excess 
feed stock is blended off or is downgraded. If the unit is oversized, 
the slack processing capacity is reported. Both situations are handled 
automatically by the program logic. 



With regard to the refinery process unit capacities, it would be 
unusual for a specified unit throughput to exactly match the volume of 
an intermediate refinery stream. If the unit is undersized, the excess 
feed stock is blended off or is downgraded. If the unit is oversized, 
the slack processing capacity is reported. Both situations are handled 
automatically by the program logic. 
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3.0 REQUIRED MODEL INPUTS 

3.1 General Rules for Fortran Namelist Input 

Punched card input is read directly by the Fuel Blending Model using 
the Fortran Namelist feature. This feature allows for a "free format" 
type of input whereby the data may be input in any sequence. However, 
certain input rules must be observed and they are listed below. 

The first card heading up the input for a new run or case 
must contain an ampersand (&) in column 2 followed by the 
name of the data list (see Section 3.2 for the instruc- 
tions for specifying data lists for the Gordian Simulation 
Model ) . 

A variable name must be completed on a card, but may begin 
in any location from column 2 onwards. Card column 1 is 
never to be used for any purpose. 

More than one variable name and more than one data entry 
may appear on a card, the only requirement being that they 
be separated by commas. 

For example: 

FCC=30000.0, GAS0IL=2000.0, GOHYC=5000.0, M0DE=2, is 
a legitimate card entry with FCC starting in card column 2 
onwards. As shown above, the = sign must immediately fol- 
low the symbolic variable name and there may be no blank 
space appearing within the symbolic variable name. However, 
the data values following the equal sign may be spaced on 
the card, as long as the last item entered on a data card 
is always a comma. 
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For example: 

FCC= 30000.0, 6AS0IL= 2000.0, 

G0HYC= 5000.0, M0DE= 2, 

are legitimate entries for two successive input cards. 
Elements of input arrays, such as the array of crude oil 
volume feed rate, must be input with the array element 
number indicated in parentheses. For example: 

CV(8)= 50000.0, CV(10)= 40000.0, 
would indicate fifty and forty thousand barrels per day 
of crude oils #8 and #10, respectively. No intervening 
blanks should appear in designating an input variable of 
this type. For example, do not input 
CV ( 8) = 50000.0, 

Data types must be respected. For example: 

M0DE=2, ICASE=1 , TITLE= 'TEST RUN', FCC= 30000.0, 

CV (1 )=50000. 0, is correct in every respect. The four 
variable types used in this program are: 

. Literal variables, which appear between quotation 
marks and which may consist of any combination of 
alphabetic or numeric information. The only variable 
of this type appearing in this program is TITLE, 
which is used to provide a run or case title. 

. Integer variables (such as ICASE and MODE) which are 
always followed by whole numbers, without a decimal 
point. 
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. Real variables (such as FCC=30000.0, VG0FCC=0.45,) 
which are followed by numbers containing decimal 
points. 

. Real variables, which are also array elements (such 
as CV(1 )=50000.0, CV(10)=20000.0, CA(2,22)=0.03 
The last card denoting the end of an input data list must 
contain &END, beginning in column 2 with no intervening 
blanks. 

3.2 Detailed Description of Specific Inputs 

The following section describes the program input using the FORTRAN 
namelist feature, item by item. If a particular input, or a group of 
inputs, is not relevant to the problem at hand, they may be omitted. A 
list of the program default values for each specific input is given in 
Section 3.3. With a few exceptions, the program default values are usu- 
ally zero. 

In this section the data type and the appropriate units are given 

for each input variable. The following abbreviations are used for units 

BPD Barrels per day. These are on a calendar 

day basis (BPCD). The program makes stream 
day adjustments based on the value of the 
input stream factor STRF. 

MMSCFD Millions of standard cubic feet per day. 

The input data stream sequence is organized as follows: 

&CDATA ... crude oil assay data base changes (optional) 

&END 

&UDATA ... process unit data base changes (optional) 

SEND 
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&UDATA1 ... process unit data base changes (optional) 

SEND 

&PDATA ... input defining the refinery feed, product and process 

SEND unit size characteristics (required) 

While crude oil and process unit data base changes are optional, 
the header cards (which contain an ampersand in column 2) are required 
at the beginning of the input data stream run sequence, even if no 
change cards are input. 

The input data items appearing between the headers &PDATA and SEND 
are described in sections 3.2.1, 3.2.2, and 3.2.7 to 3.2.9. Default 
values are automatically established within the program, as discussed in 
section 3.3; hence input data are required only where different from the' 
default values. The optional crude oil assay data items appearing 
between the headers &CDATA and SEND are described in section 3.2.3. 
Optional process data change items appearing between &UDATA and SEND are 
those described in sections 3.2.4 and in section 3.2.5 except for 
viscosity and research octane variables. Optional process data change 
items appearing between &UDATA1 and SEND are the viscosity and octane 
variables and the items described in section 3.2.6. 
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3.2.1 

Title and Case Study Feature 


Name 

Type Uni ts 

Description 


TITLE 

Literal 

The run or case title which will appear at the 
top of each output page. Up to 70 alpha-numeric 
characters may be used. 

ICASE 

Integer 

The case study indicator, 
assumed that the run is a 

If not entered it is 
"stand alone run" with 


no fixed relation to preceding or succeeding runs 
Input ICASE=1 for the starting base case of a 
case study sequence. 

Input ICASE=2 for a subject case to be dependent 
on the starting case. 

Input ICASE=3 for the final case of a case study 
sequence. 

A case study sequence may be followed by stand 
alone runs or a new case study sequence can be 
initiated, if desired. 

The program logic is such that all subject case 
changes are made to the starting base case of a 
case study sequence. Therefore, subject case 
changes are not additive in nature. 


3.2.2 Process Unit Capacities, Modes of Operation 
and Internal Stream Disposition ■ 


Name 

1m. 

Uni ts 

Description 

CV(N) 

Real 

Array 

BPD 

The crude unit charge rate of crude N, where 
the code integer N applies to a specific crude. 
A list of those crudes currently represented 
in the model is given in Exhibit 4, attached. 

VACD 

Real 

BPD 

The vacuum distillation capacity. For petro- 
leum-based, expressed as reduced crude feedrate. 

VACS 

Real 

BPD 

The vacuum distillation capacity for shale 
oil based crude. 

VACC 

Real 

BPD 

The vacuum distillation capacity for coal oil 
based synthetic crude. 

FCC* 

Real 

BPD 

The refinery catalytic cracking capacity, ex- 
pressed as total gas oil feedrate (fresh feed). 

CONV 

Real 

Volume 

Percent 

The conversion (percentage gas oil disappearance) 
of the fluid catalytic cracking unit. 


* Presently available data is limited to petroleum based fractions. The shale and coal 
oil h.ydroprocessing_units should be specified for process ing_those i species. 
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Name 

Jm. 

Uni ts 

Description 

ALKY 

Real 

BPD 

The refinery alkylation unit capacity expressed 
as total alkylate product. 

POLY 

Real 

BPD 

The refinery polymerization unit capacity 
expressed as polymer product rate. 

MDDES* 

Real 

BPD 

The refinery middle distillate desulfurizer 
feed capacity expressed as total feed rate. 

GODES* 

Real 

BPD 

The refinery gas oil desulfurizer feed capacity 
expressed as total feed rate. 

SOKHT 

Real 

BPD 

The refinery shale oil kerosene hydrotreater 
capacity expressed as kerosene feed rate. 

SVSKHT 

Integer 

— 

The designated mode of operation of the shale 
oil kerosene hydrotreater. Enter 1 for mild, 

2 for intermediate, and 3 for severe operation. 

COKHT 

Real 

BPD 

The refinery coal oil kerosene hydrotreater 
capacity expressed as kerosene feed rate. 

SVCKHT 

Integer 

— 

The designated mode of operation of the coal 
oil kerosene hydrotreater. Enter 1 for mild, 

2 for intermediate, and 3 for severe operation. 

SGOHYC 

Real 

BPD 

The refinery shale oil gas oil hydrocracker 
capacity expressed as total feed rate. 

SVSGHT 

Integer 

*• 

The designated mode of operation of the shale 
oil gas oil hydrocracker. Enter 1 for mild, 

2 for intermediate, and 3 for severe operation. 

CGOHYC 

Real 

BPD 

The refinery coal oil gas oil hydrocracker 
capacity expressed as total feed rate. 

SVCGHT 

Integer 

— 

The designated mode of operation of the coal 
oil gas oil hydrocracker. Enter 1 for mild, 

2 for intermediate, and 3 for severe operation. 

COKER 

Real 

BPD 

The refinery coker unit capacity expressed as 
vacuum bottoms feed rate. 

VB 

Real 

BPD 

The refinery visbreaker unit capacity expressed 
as vacuum bottoms feed rate. 

TCC 

Real 

BPD 

The refinery thermal cracker unit capacity 
expressed as gas oil feed rate. 


Specify these units for processing petroleum based fractions. Specify 
the shale and coal oil hydroprocessing units for processing those species. 
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Name 

Jm. 

Units 

Description 

MDHYC* 

Real 

BPD 

The refinery middle distillate hydrocracker 
capacity expressed as total feed rate. 

60HYC* 

Real 

BPD 

The refinery petroleum gas oil hydrocracker capacity 
expressed as total feed rate. 

MODE 

Integer 


The designated mode of operation of the refinery 
gas oil hydrocracker. Enter M0DE=1 for a 
gasoline operatidn and M0DE=2 for a middle 
distillate operation. 

BI 

Real 

BPD 

The refinery butane isomerization capacity 
expressed as the normal butane feed rate. 

HYD 

Real 

MMSCFD 

The refinery hydrogen unit capacity expressed 
as the hydrogen production rate. 

REF 

Real 

BPD 

The refinery catalytic naphtha reformer capacity 
expressed as total naphtha feed rate. 

REFRON 

Real 

Research 

Octane 

No., 

Clear 

The catalytic reformer severity expressed as 
the clear research octane number of the total 
platformate. 

HTKERO 

Real 

BPD 

The refinery petroleum kerosene hydrotreater capacity 


expressed as the kerosene feed rate. 


Internal Stream Dispositions 

If alternative destinations exist for any of the internal refinery 
streams, the following inputs serve to define the fraction of the stream 
being routed to each alternative process unit (values for these variables 
should not be input for a given stream if the refinery configuration is 
such that multiple destinations do not exist). The programmed default 
is to prorate stream volume by the ratio of available downstream process 
unit capacities. 


* Specify these units only for processing petroleum based fractions. Specify 
the shale and coal oil hydroprocessing units for processing those species. 
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If the total of the sum of the volume fractions going to alternative 
destinations for a particular stream are specified at less than 1.0, the 
balance of that stream is assumed to go directly to fuel oil blending. 

For the low-boiling internal streams (methane, ethane, propane, 
butanes, and the corresponding olefins), distribution priorities are 
prescribed. Alkylation has the first priority for isobutane and butylene. 
The polymerization unit uses propylene first, and if additional capacity 
is left, is uses the excess butylene. For hydrogen production, propane 
is used first and then ethane and lighter. The overall hydrogen require- 
ments of the refinery are filled by first the hydrogen produced by the 
catalytic reformer and then, if needed, by hydrogen from the hydrogen 
plant. The balance of the low-boiling streams not used internally goes 
to fuel gas and LPG. 


Name 

Jm. 

Units 

Description 

DSLDES 

Real 

Vol ume 
Fraction 

The volume fraction of the refinery diesel 
stream going to the middle distillate 
desulfurizer. 

DSLMDH 

Real 

Vol ume 
Fraction 

The volume fraction of the refinery diesel 
stream going to the middle distillate hydro- 
cracker. 

LCODES 

Real 

Volume 

Fraction 

The volume fraction of light catalytic cycle 
oil going to the middle distillate desulfurizer. 

LCOMDH 

Real 

Volume 
Fracti on 

The volume fraction of light catalytic cycle 
oil going to the middle distillate hydrocracker. 

VGODES 

Real 

Volume 

Fraction 

The volume fraction of the refinery gas oil 
stream going to gas oil desulfurization. 

VGOHYC 

Real 

Vol ume 
Fraction 

The fraction of the above stream going to gas 
oil hydrocracking. 

VGOFCC 

Real 

Vol ume 
Fraction 

The fraction of the above stream going to gas 
oil catalytic cracking. 

VGOTCC 

Real 

Vol ume 
Fraction 

The fraction of the above stream going to gas 
oil thermal cracking. 
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CGODES 

COGOHY 

CGOFCC 

Real 

Volume 

Fractions 

These three input variables denote the volume 
fraction of total coker gas oil to the respec- 
tive destinations listed above. 

VBTMVB 

Real 

Vol ume 
Fraction 

The volume fraction of vacuum bottoms going to 
the visbreaker. 

VBTMCO 

Real 

Volume 

Fraction 

The volume fraction of vacuum bottoms going to 
the coker. 

SLKKHT 

• 

Real 

Volume 

Fraction 

The volume fraction of light kerosene (400 to 
525°F EP) going to the shale oil kerosene 
hydrotreater. 

SHKKHT 

Real 

Volume 

Fraction 

The volume fraction of heavy kerosene (525 to 
650°F EP) going to the shale oil kerosene 
hydrotreater. 

CLKKHT 

Real 

Volume 

Fraction 

The volume fraction of light kerosene (400 to 
525° F EP) going to the coal oil kerosene 
hydrotreater. 

CHKKHT 

Real 

Volume 

Fraction 

The volume fraction of heavy kerosene (525 
to 650°F EP) going to the coal oil kerosene 
hydrotreater. 

SGOHC 

Real 

Volume 

Fraction 

The volume fraction of gas oil going to the 
shale oil gas oil hydrocracker. 

CGOHC 

Real 

Volume 

Fraction 

The volume fraction of gas oil going to the 
coal oil gas oil hydrocracker. 

3.2.3 

Crude Assay Information 

The 

following 

input variables relate to the crude oil assay information 


which may be entered directly as input. The values may apply to entirely 
new feeds or may be used to override crude oil assay information currently 


stored in the model. The crude assay data changes are entered first in 
the sequence between the &CDATA and SEND header cards. 

Crude oil assay code numbers 1 through 26 are reserved for crude oil 
assays presently stored within the program. Code numbers 27 through 35 may 
be used to describe crude oils and feeds which are not in the permanent 
data base. A list of data base crude oils is given in Exhibit 4. 
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Refinery charge streams such as gas oils may be specified by entering 
a fractional gas oil yield of 1.0 for the crude assay, along with the gas 
oil properties. The remaining crude oil assay values are not to be input 
in this instance since the gas oil data has fully described the refinery 
charge stream. This input data philosophy may be extended to any type of 
refinery charge stream, naphtha, kerosene, diesel, etc., which is imported 
to the refinery for processing within the normal plant sequence. 


Name 

Im. 

Units 

Description 

CA(N,1) 

Real 

Array 

Vol ume 

Fraction 

(BBL/BBL) 

Crude assay yield of ethane and lighter for 
crude N. 

C A ( N , 2 ) 

II 

(I 

Yield of propane. 

CA(N,3) 

II 

II 

Yield of isobutane. 

CA(N,4) 

II 

II 

Yield of normal butane. 

CA(N,5) 

II 

II 

Yield of pentane. 

CA(N,6) 

It 

II 

Yield of light straight run gasoline. 
(IBP-250°F) 

CA(N,7) 

II 

11 

Yield of heavy straight run naphtha- 
reformer feed. (250 to 400° F) 

CA(N,8) 

II 

II 

Yield of light kerosene (400 to 525°F). 

CA(N,9) 

II 

II 

Yield of heavy kerosene (525 to 650°F). 

C A ( N , 1 0 ) 

11 

II 

Yield of vacuum gas oil (650 to 1,050°F EP). 

CA(N ,11 ) 

11 

II 

Yield of vacuum bottoms (1,050°F +). 

C A ( N , 1 2 ) 

II 

Degrees 

API 

Gravity 

Crude oil API gravity. 

1 

CA (N , 13) 

II 

ll 

Light straight run API gravity. 

CA(N,14) 

II 

11 

Heavy straight run API gravity. 

CA(N,15) 

11 

11 

Light kerosene API gravity. 

CA(N , 16) 

II 

ll 

Heavy kerosene API gravity. 
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CA(N,17) 

CA(N,18) 

CA(N,19) 

CA(N,20) 

CA( N , 21 ) 

CA(N,22) Real 
Array 

CA(N,23) 

CA(N,24) 

CA(N,25) 

CA(N,26) 

CA(N,27) 

CA(N,28) 

CA(N,29) 

CA(N,30) 

CA(N, 31 ) 
CA(N,32) 
CA(N,33) 
CA(N,34) 
CA(N,35) 
CA(N,36) 
CA(N,37) 
CA(N,38) 
CA(N,39) 
CA(N,40) 

CA(N,41 ) 
CA(N,42) 


" Vacuum gas oil API gravity. 

" Vacuum bottoms API gravity. 

Weight Crude oil sulfur content. 

Percent 

" Light straight run sulfur content. 

" Heavy straight run sulfur content. 

Weight Light kerosene sulfur content. 

Percent 

" Heavy kerosene sulfur content. 

" Vacuum gas oil sulfur content. 

" Vacuum bottoms sulfur content. 

" Crude oil nitrogen content. 

" Light straight run nitrogen content. 

" Heavy straight run nitrogen content. 

" Light kerosene nitrogen content. 

" Heavy kerosene nitrogen content. 

" Vacuum gas oil nitrogen content. 

" Vacuum bottoms nitrogen content. 

" Crude oil hydrogen content. 

" Light straight run hydrogen content. 

" Heavy straight run hydrogen content. 

" Light kerosene hydrogen content. 

" Heavy kerosene hydrogen content. 

" Gas oil hydrogen content. 

" Vacuum bottoms hydrogen content. 

Viscosity Light kerosene viscosity index at 210°F. 

Blending 

Index 


Heavy kerosene viscosity index at.210°F 

Gas oil viscosity index at 210°F. 
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II 


CA(N,43) 

11 

11 

Vacuum bottoms 

viscosity index at 210°F. 

CA(N,44) 

Real 

Array 

Weight 

Percent 

Light 

kerosene 

paraffin content. 

CA(N,45) 

ll 

ll 

Light 

kerosene 

naphthene content. 

CA(N,46) 

li 

ll 

Light 

kerosene 

aromatics content. 

CA(N,47) 

ll 

il 

Heavy 

kerosene 

paraffin content. 

CA(N,48) 

it 

ll 

Heavy 

kerosene 

naphthene content. 

CA(N,49) 

ll 

ll 

Heavy 

kerosene 

aromatics content. 

CA(N,50) 

li 

?F 

Light 

kerosene 

freezing point. 

CA(N,51 ) 

ll 

°F 

Heavy 

kerosene 

freezing point. 

CA(N,52) 

ll 

mm 

Light 

kerosene 

smoke point. 

CA(N,53) 

ll 

mm 

Heavy 

kerosene 

smoke point. 

CA(N,54) 

ll 

BTU/# 

Crude 

oil heat 

of combustion. 

CA(N,55) 

ll 

H 

Light 

kerosene 

heat of combustion. 

CA(N,56) 

ll 

it 

Heavy 

kerosene 

heat of combustion. 

CA(N,57) 

ll 

* 

Species indicator 

1.0 = Coal Oil 

2.0 = Shale Oil 

3.0 = Petroleum 


Note : N corresponds to the crude oil code number given in Exhibit 

4. The boiling point distillation ranges of the liquid 
crude oil cuts are given in parentheses and correspond to 
ASTM D-86 boiling points. The Viscosity Blend Index 
referred to is described in Exhibit 6. 

3.2.4 Process Unit Data 

The following input variables affect the internally stored program 
values of the process unit yields and utility consumptions. It is only 
necessary to input these variables when it is desired to replace the 
data stored for any of the process units with values of the user's 
choice. Changes to the process unit data base are entered second in the 
input sequence, after the crude assay data base changes between &UDATA 
and SEND header cards. 
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All of the input variables listed below are elements of real arrays. 
The units are fractional barrels of the particular stream yielded per 
barrel of the indicated feed stream. Coke, hydrogen, refinery gas, and 
process unit steam and fuel consumption are in barrels of fuel oil 
equivalent per barrel of feed.* Electricity usage is in kilowatt hours 
(KWH) per barrel of process unit feed. If the process unit produces 
steam as opposed to using it, a negative number should be entered for 
steam consumption. 

As may be seen from analyzing the input variable names listed below, 
all of the product streams have been assigned characteristic code numbers. 
These are listed in Exhibit 5. 

Crude Unit, 

Per Barrel of Crude 


Oil Feed 

Petroleum 

Shale Oil 

Coal Oil 

Steam consumption 

CUPYLD(l) 

CUSYLD(l) 

CUCYLD(l) 

Fuel consumption 

CUPYLD(2) 

CUSYLD(2) 

CUCYLD(2 

Electricity consumption 

CUPYLD(3) 

CUSYLD(3) 

CUCYLD(3) 


Visbreaker 


Coker 


per barrel 

of vacuum bottom 

feed 

Steam consumption 

VBYLD(l) 


COYLD(l) 

Fuel consumption 

VBYLD(2) 


C0YLD(2) 

Electricity .consumption 

VBYLD(3) 


C0YLD(3) 

Refinery gas yield 

VBYLD(5) 


C0YLD(5) 

Propane yield 

. VBYLD 6 


C0YLD(6) 

Propylene yield 

VBYLD(7) 


C0YLD(7) 

Isobutane yield 

VBYLD(8) 


C0YLD(8 

Normal butane yield 

VBYLD 9) 


C0YLD(9) 

Butylene yield 

VBYLD(IO) 


COYLD(IO) 

Light gasoline 

VBYLD (31) 


C0YLD( 31 ) 

Heavy gasoline 

VBYLD( 32) 


C0YLD(32) 

Total gas oil 

VBYLD (35) 


C0YLD(35) 

Coke 

- 


C0YLD(36) 

Pitch 

VBYLD (51) 


- 


6,050,000 BTU per barrel fuel oil equivalent (BF0E), corresponding to 
the net heating value. 
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Catalytic Cracker Yields 

@ low conversion 
per barrel of 

@ high conversion 
vacuum gas oil feed 

Steam consumption 

CGYLDL(l) 


CGYLDH (1 ) 

Fuel consumption 

CGYLDL(2) 


CGYLDH (2 

Electricity consumption 

CGYLDL(3) 


CGYLDH (3) 

Refinery gas yield 

CGYLDL(5) 


CGYLDH(5) 

Propane yield 

CGYLDL(6) 


CGYLDH(6) 

Propylene yield 

CGYLDL(7) 


CGYLDH(7) 

Isobutane yield 

CGYLDL(8) 


CGYLDH(8) 

Normal butane yield 

CGYLDL(9) 


CGYLDH(9) 

Butylene yield 

CGYLDL(IO) 


CGYLDH (10) 

Cracked gasoline 

CGYLDL(37) 


CGYLDH(37) 

Light cycle oil 

CGYLDL(38) 


CGYLDH(38) 

Cracker bottoms 

CGYLDL(39) 


CGYLDH(39) 


per barrel 

of coker gas oil feed 

Steam consumption 

CCYLDL(l) 


CCYLDH ( 1 ) 

Fuel consumption 

CCYLDL(2) 


CCYLDH(2) 

Electricity consumption 

CCYLDL (3) 


CCYLDH(3) 

Refinery gas yield 

CCYLDL(5) 


CCYLDL(5) 

Propane yield 

CCYLDL (6) 


CCYLDH(6) 

Propylene yield 

CCYLDL(7) 


CCYLDH (7) 

Isobutane yield 

CCYLDL (8) 


CCYLDH(8) 

Normal butane yield 

CCYLDL (9) 


CCYLDH(9) 

Butylene yield 

CCYLDL (10) 


CCYLDH (10) 

Cracked gasoline 

CCYLDL (37) 


CCYLDH(37) 

Light cycle oil 

CCYLDL (38) 


CCYLDH (38) 

Cracker bottoms 

CCYLDL (39) 


CCYLDH(39) 


CONVL Low Catalytic Cracker conversion — volume 
percent of gas oil disappearance. 

CONVH High Catalytic Cracker conversion — volume 
percent of gas oil disappearance 


Thermal Cracker Yields 


per barrel of vacuum gas oil feed 


Steam consumption 

TCCYLD(l) 

Fuel Consumption 

TCCYLD(2) 

Electricity consumption 

TCCYLD(3 

Refinery gas yield 

TCCYLD(5) 

Propane yield 

TCCYLD(6) 

Propylene yield 

TCCYLD(7) 

Isobutane yield 

TCCYLD(8) 

Normal Butane yield 

TCCYLD(9) 

Butylene yield 

TCCYLD(IO) 

Cracked gasoline 

TCCYLD(37) 

Thermal fuel 

TCCYLD( 52 ) 
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Alkylation Unit Yields 

per barrel of butylene feed 


Steam consumption 
Fuel consumption 
Electricity consumption 
Iso Butane consumption 
Light Alkylate yield 
Heavy Alkylate yield 


ALYLD(l) 

ALYLD(2) 

ALYLD(3) 

ALYLD(8) (Enter a negative number) 
ALYLD( 21 ) 

ALYLD(22) 


Polymerization Unit Yields 


per barrel of 
butylene feed 


per barrel of 
propylene feed 


Steam consumption POYLDB(l) 
Fuel consumption P0YLDB(2) 
Electricity consumption P0YLDB(3) 
Normal butane yield P0YLDB(9) 
Polymer gasoline yield P0YLDB(23) 


POYLDP(l) 
POYLDP ( 2 ) 
POYLDP ( 3 ) 
POYLDP ( 9 ) 
POYLDP (23) 


Middle Distillate Desulfurizer Yields 



per barrel of 

per barrel of light 


diesel feed 

. cycle oil feed 

Steam consumption 

DDYLDD(l) 

DDYLDL(l) 

Fuel consumption 

DDYLDD(2) 

DDYLDL(2) 

Electricity consumption 

DDYLDD(3) 

DDYLDL(3) 

Hydrogen consumption* 

DDYLDD(4) 

DDYLDL(4) 

Refinery gas yield 

DDYLDD(5) 

DDYLDL(5) 

Desulfurized product yield 

DDYLDD 41) 

DDYLDL 40) 

Gasoline yield 

DDYLDD(45) 

DDYLDL(45) 

Gas Oil 

Desulfurizer Yields 



per barrel of 

per barrel of total 


vacuum gas oil 

coker gas oil 


Steam consumption GDYLDG(l) GDYLDC(l) 
Fuel consumption GDYLDG(2) GDYLDC(2) 
Electricity consumption GDYLDG(3) GDYLDC(3) 
Hydrogen consumption* GDYLDG(4) GDYLDC(4) 
Refinery gas yield GDYLDG(5) GDYLDC(5) 
Distillate product yield GDYLDG(43) GDYLDC(46) 
Desulfurized gas oil yield GDYLDG(44) GDYLDC(47) 


Enter a negative number for hydrogen in all processes where 
hydrogen is consumed. 
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Middle Distillate Hydrocracker Yields 


Steam consumption 
Fuel .consumption 
Electricity consumption 
Hydrogen consumption* 
Refinery gas yield 
Propane yield 
Iso Butane yield 
Normal Butane yield 
Pentanes yield 
Hydrocrackate yield 


per barrel of 
diesel . feed 

per barrel of light 
cycle oil feed 

DHYLDD(l) 

DHYLDL(l) 

DHYLDD(2) 

DHYLDL(2) 

DHYLDD(3) 

DHYLDL(3) 

DHYLDD(4) 

DHYLDL(4 

DHYLDD(5) 

DHYLDL(5) 

DHYLDD(6) • 

DHYLDL(6) 

DHYLDD(8) 

DHYLDL(8) 

DHYLDD(9) 

DHYLDL(9) 

DHYLDD(13) 

DHYLDL( 13) 

DHYLDD(14) 

DHYLDL(14) 


Gas Oil Hydrocracker Yields 


Steam consumption 
Fuel consumption 
Electricity consumption 
Hydrogen consumption* 
Refinery fuel gas yield 
Propane yield 
Iso Butane yield 
Normal Butane yield 
Light Hydrocrackate yield 
Heavy Hydrocrackate yield 
Distillate yield 


per barrel of vacuum 

gas oil feed 

Gasoline Distillate 
Mode Mode 


GGYLDG(l) GMYLDG(l) 
GGYLDG(2) GMYLDG(2) 
GGYLDG(3) GMYLDG(3) 
GGYLDG(4) GMYLDG(4) 
GGYLDG(5) GMYLDG(5) 
GGYLDG(6) GMYLDG(6) 
GGYLDG(8) GMYLDG(8) 
GGYLDG(9) GMYLDG(9) 
GGYLDG(15) GMYLDG(15) 
GGYLDG(14) GMYLDG(14) 
GGYLDG(17) GMYLDG(17) 


per barrel of total 
coker gas oil feed 
Gasoline Distillate 
Mode Mode 


GGYLDC(l) GMYLDC(l) 
GGYLDC(2) GMYLDC(2) 
GGYLDC(3) GMYLDC(3) 
GGYLDC(4 GMYLDC(4) 
GGYLDC ( 5 ) GMYLDC(5) 
GGYLDC(6) GMYLDC(6) 
GGYLDC(8) GMYLDC(8) 
GGYLDC(9) GMYLDC(9) 
GGYLDC (15) GMYLDC(15) 
GGYLDC( 14) GMYLDC (1 4) 
GGYLDC (48) GMYLDC (48) 


Enter a negative number for hydrogen in all processes where hydrogen 
is consumed. 


33 


Catalytic Reformer Yields 


Steam consumption 
Fuel consumption 
Electricity consumption 
Hydrogen yield* 

Refinery gas yield 
Propane yield 
Iso Butane yield 
Normal Butane yield 
Light Platformate yield 
Heavy Platformate yield 
C-5, C-6 Feed Prep yield 
Heavy-Naptha (Feed 
Prep yield) 


per barrel of naththa feed 


0 Low Severity 

RYLDLN(l) 
RYLDLN(2) 
RYLDLN(3) 
RYLDLN (4) 
RYLDLN(5) 
RYLDLN (6) 
RYLDLN(8) 
RYLDLN (9) 
RYLDLN(18) 
RYLDLN (19) 
RYLDLN (12) 

RYLDLN(16) 


0 High Severity 

RYLDHN(l) 

RYLDHN(2) 

RYLDHN (3) 

RYLDHN 4) 
RYLDHN(5) 

RYLDHN (6) 

RYLDHN (8) 

RYLDHN (9) 
RYLDHN(18) 
RYLDHN (19) * 
RYLDHN (12) 

RYLDHN( 1 6) 


For Heavy hydrocrackate reformer feed, use the same input sequence but with 
RYLDLH( J and RYLDHH(J. 


For Heavy coker or visbreaker naphtha reformer feed, use the same input 
sequence but with RYLDLC(_) and RYLDHC(_). 


RONLOW The debutanized total platformate clear research octane 
number corresponding to low reformer severity. 

RONHI The debutanized total platformate clear research octane 
number corresponding to high reformer severity. 


Hydrogen Plant Yields 


per barrel FOE per liquid bbl 

from refinery gas from propane 


Steam consumption 
Fuel consumption 
Electricity consumption 
Hydrogen production** 


HYYLDG(l) 
HYYLDG(2) 
HYYLDG(3) 
HYYLDG (4) 


HYYLDP(l) 

HYYLDP(2) 

HYYLDP(3) 

HYYLDP(4) 


i 


* Enter a positive number here in barrels FOE hydrogen per barrel of 
feed stream (95% pure hydrogen) since there is a production rather 
than a consumption of hydrogen. 

** A positive number. 
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Butane Isomerizer 


per barrel normal butane feed 


Steam consumption 

BIYLD(l) 

Fuel consumption 

BIYLD(2) 

Electricity consumption 

BIYLD(3) 

Refinery gas yield 

BIYLD(5) 

Propane yield 

BIYLD(6) 

Iso Butane yield 

BIYLD(8) 

Pentanes yield 

BIYLD(ll) 


Petroleum Kerosene Hydrotreater Yields 



per barrel of kerosene feed 

Steam consumption 

HKYLD(l) 

Fuel consumption 

HKYLD(2) 

Electricity consumption 

HKYLD(3) 

Hydrogen consumption 

HKYLD(4) 

Refinery gas yield 

HKYLD(5) 

Hydrotreater gasoline yield 

HKYLD(45) 

Hydrotreater kerosene 

HKYLD(50) 


Shale Oil Kerosene Hydrotreater Yields 


per bbl of light per bbl of heavy 

kerosene feed kerosene feed 

(400 to 525 F EP) (525 to 650°F EP ) 


Low Severity 


Low Severity 


Steam consumption SSKYLl(l)* 

Fuel consumption SSKYL1 (2 ) 

Electricity consumption SSKYLl(3) 

Hydrogen consumption* SSKYL1(4) 

Refinery fuel gas yield SSKYL1(5) 

Hydrotreated kerosene SSKYL1(58) 


SSKYH1 (1 )* 
SSKYH1 (2) 
SSKYH1 (3) 
SSKYH1 (4) 
SSKYH1 (5) 
SSKYH1 (59) 


Note: For intermediate and high severity operations, replace 

the.. numeral 1 (i .e. , the sixth character) by 2 Or 3, 
respectively. 
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Coal Oil Kerosene Hydrotreater Yields 



per bbl of light 

per bbl of heavy 


kerosene feed 

kerosene feed 


(400 to 525°F EP) 

(525 to 650°F EP) 


Low Severity 

Low Severity 

Steam consumption 

C0KYL1 (1 )* 

C0KYH1 (1 )* 

Fuel consumption 

C0KYL1 (2) 

C0KYH1 (2) 

Electricity consumption 

C0KYL1 3 

C0KYH1 3 

Hydrogen consumption** 

C0KYL1 (4) 

C0KYH1 (4) 

Refinery fuel gas yield 

C0KYL1 (5) 

C0KYH1 (5) 

Hydrotreated kerosene 

C0KYL1 (64) 

C0KYH1 (65) 

Shale and Coal 

Oil for Oil Hydrocrackers 



Shale Oil 

Coal Oil 


Low Severity 

Low Severity 

Steam consumption 

SSGY1 (1 )* 

C0GY1 (1 )* 

Fuel consumption 

SSGY1 (2) 

C0GY1 (2) 

Electricity consumption 

SSGY1 (3) 

C0GY1 (3) 

Hydrogen consumption* 

SSGY1 (4) 

C0GY1 (4) 

Refinery fuel gas yield 

SSGY1 (5) 

C0GY1 (5) 

Propane yield 

SSGY1 (6) 

C0GY1 (6) 

Isobutane yield 

SSGY1 (8) 

C0GY1 (8) 

N- butane yield 

SSGY1 (9) 

C0GY1 (9) 

Light hydrocrackate 



(250°F EP) 

SSGY1 (15) 

C0GY1 (15) 

Heavy hydrocrackate 



(400°F EP) 

SSGY1 (14) 

C0GY1 (14) 

Light hydrotreated kerosene 



(400° - 525°F EP) 

SSGY1 (58) 

C0GY1 (64) 

Heavy hydrotreated kerosene 



(525° - 650°F EP) 

SSGY1 (59) 

C0GY1 (65) 

Hydrotreated gas oil 



(650°F +) 

SSGY1 (60) 

C0GY1 (66) 


* Note : For intermediate and high severity operations, replace 

the numeral 1 ( i .e . , the fifth character) by 2 or 3, 
respectively. 

3.2.5 Intermediate Stream Properties 

The following input variables relate to the properties of the indi- 
vidual streams which are yielded from the various refinery process 
units. They are not required as input unless it is desired to replace 
the data base values contained within the program. Input items described 


** Enter a negative number. 
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in this section are included between the headers &UDATA and SEND, except 
for viscosity (VIS or V) and octane (OC) properties, which are included 
between &UDATA1 and &END. 

For aviation turbine fuel and fuel oil blending components the pro- 
perties include the API gravity, sulfur content, nitrogen content and the 
viscosity blending index value. Hydrogen contents are calculated 
from a correlation of specific gravity and mid-boiling points except for 
select streams where a hydrogen content may be specified to override the 
calculated value. Gasoline blending component properties also include 
the Research (blending) Octane Number for 0 and 3 ml of TEL. 

The viscosity blending numbers used must correspond to Saybolt 
Seconds at 210 degrees Fahrenheit so that they are consistent with the 
values contained in the model. A viscosity conversion chart is attached 
as Exhibit 6. 

With the exceptions to be noted, the input variables are all elements 
of real arrays. The butane blending octane numbers are (non-array) real 
numbers. 

A consistent logic applies to the naming of intermediate stream 


properties. This logic is 

illustrated 

below using the kerosene hydro- 

treater unit as an example: 



Description 

Array 

Remarks 

Yield Array 

HKYLD(I) 

See section 2.2.4 for details. 

API Gravity 

HKAPI(I) 

For pentanes and heavier' 
streams. 

Sulfur Distribution 

HKSUL(I) 

Fractional sulfur distribution 
among products.* 

Nitrogen Distribution 

HKNIT(I) 

Fractional nitrogen distributi 
among products.* 


These should add to 1.0, unless there are external losses, such as 
loss of sulfur from the combustion of catalytic coke. 
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Hydrogen Content HKHYD(I) Hydrogen content, wt. % 

Viscosity Blending Number HKVIS(I) For kerosene and heavier streams, 

(input between &UDATA1 and &END) 

The above logic is used to name the properties of streams for 
all of the process units using the three-character set YLD and which are 
ennumerated in Section 3.2.4. The character set YLD is replaced by API, 

SUL, NIT, HYD or VIS to input the appropriate property. The letter "I" 
shown under the heading Array refers to the stream identification number 
provided in Section 3.2.4 for process units employing a single character 
to describe the yield array. The letter Y is replaced by A, S, N and H 
to input API, sulfur distribution, nitrogen distribution and hydrogen 
content, respectively, for units employing a single character. These 


are listed below: 


Uni t 

Array 


Shale Oil Kerosene Hydrotreater 
Low Severity - Light Kerosene 
Medium Severity - Light Kerosene 
High Severity - Light Kerosene 
Low Severity - Heavy Kerosene 
Medium Severity - Heavy Kerosene 
High Severity - Heavy Kerosene 

SSKYL1 (I) 
SSKYL2 ( I ) 
SSKYL3(I j 
SSKYH1 ( I ) 
SSKYH2(I] 
SSKYH3(I) 


Coal Oil Kerosene Hydrotreater 
Low Severity - Light Kerosene 
Medium Severity - Light Kerosene 
High Severity - Light Kerosene 
Low Severity - Heavy Kerosene 
Medium Severity - Heavy Kerosene 
High Severity - Heavy Kerosene 

C0KYL1 ( I ] 
C0KYL2(i; 
C0KYL3(i; 

cokyhki; 

C0KYH2(I1 

C0KYH3(i; 

1 

1 

1 

) 

I 

) 

Shale Oil Gas Oil Hydrocracker 
Low Severity 
Medium Severity 
High Severity 

SSGYl(I) 

SSGY2(I) 

SSGY3(I) 

i 

Coal Oil Gas Oil Hydrocracker 
Low Severity 
Medium Severity 
High Severity 

COGYl(I) 

C0GY2(I) 

C0GY3(I) 


For those processing units, viscosity is expressed as 

non-array 


real numbers. 
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Unit 


Variable 


Shale Oil Kerosene Hydrotreater 


Light Kerosene SSKVL1 

Heavy Kerosene SSKVH1 

Shale Oil Gas Oil Hydrocracker 

Light Kerosene SSGVL1 

Heavy Kerosene SSGVH1 

Fuel Oil SSGVG1 


For the coal oil viscosities corresponding to the above, replace SS 
with CO. 


Hydrogen content of streams may be input only for certain 
process units — generally those units which yield some portion of their 
product directly into aviation turbine fuel blends. If values are not 


input, the hydrogen content is calculated from a generalized correlation 
employing specific gravity and mid-boiling point. For units not listed 


below, hydrogen content may not be specified and the correlation results 
may not be overriden by input specification. 


List of Units for Which Product 
Stream Hydrogen Contents May be Specified 


Unit 


Array Name 


Distillate desulfurizer (petroleum) 

Straight run stock DDHYDD(I) 
FCC light cycle oil DDHYDL(I) 
Light coker/ visbreaker gas oil DDHYDC(I) 

Gas Oil hydrocracker (petroleum) 

Straight run stock-gasoline mode GGHYDG(I) 
Coker/visbreaker stock-gasoline mode GGHYDC(I) 
Straight run stock-distillate mode GMHYDG(I) 
Coker/ visbreaker stock-distillate mode GMHYDC(I) 
Kerosene hydrotreater (petroleum) HKHYD(I) 


Low Medium High 


Unit 

Severity 

Severi ty 

Severity 

Shale oil light kerosene 
hydrotreater 

SSKHLl(I) 

SSKHL2 ( I ) 

SSKHL3(I) 

Shale oil heavy kerosene 
hydrotreater 

SSKHH1 (I) 

SSKHH2 ( I ) 

SSKHH3 ( I ) 
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Coal oil light kerosene 
hydrotreater 

Coal oil heavy kerosene 
hydrotreater 

Shale oil gas oil hydrocracker 
Coal oil gas oil hydrocracker 


C0KHL1 (I) 

C0KHH1 (I) 

SSGHl(I) 

COGHl(I) 


C0KHL2(I) C0KHL3 ( I ) 

C0KHH2 ( I ) C0KHH3 ( I ) 
SSGH2(I) SSGH3 (I) 
C0GH2(I) C0GH3 (I) 


Input variable specification for gasoline blending component 


streams are given below (note that butane values are non-array). All 
these quantities are input in between &UDATA1 and SEND. 


Research Octane Number 


0 0 ml . TEL 0 3 ml . TEL 


Butanes 

Reformer Feed Prep C5/C6 

Mid. Dist. Hydrocracker C5's 

Gas Oil Hydrocracker-Light Hydrocrackate 

Crude Unit Straight Run Gasoline (incl. C5's) 

Cracker Gasoline (Full range) 

From vacuum gas oil at low conversion* 

From vacuum gas oil at high conversion 
From coker gas oil at low conversion 
From coker gas oil at high conversion 
Visbreaker/Coker light gasoline** 

Light Alkylate 
Heavy Alkylate 
Polymer gasoline 
Light platformate-low severity 
Light platformate-high severity 
Heavy platformate-high severity 
Heavy platformate-high severity 
Desulfurizer Overhead Gasoline 


C40CL C40CH 
R0CLD(12) ROCHD (12) 
DH0CL(13) DH0CH(13) 

GGOCL (15) GG0CH(15) 

CU0CL(24) CUOCH (24 ) 

CG0CLL(37) CG0CHL(37) 
CG0CLH(37) CG0CHH(37) 
CC0CLL(37) CC0CHL(37) 
CC0CLH(37) CCOCHH ( 37 ) 
COOCLL ( 31 ) C00CHL(31 ) 
ALOCL (21 ) ALOCH ( 21 ) 
AL0CL(22) ALOCH (22) 
P00CL(23) POOCH (23) 
R0CLLL(18) R0CHLL(18) 
R0CLHL(18) R0CHHL(18) 
R0CLLH(19) R0CHLH(19) 
R0CLHH(19) R0CHHH(19) 
DD0CL(45) DD0CH(45) 


* Uses the low conversion input variables for a thermal gas oil 
cracker. 

** For 250° F EP light gasoline (below reforming range) produced 
from vacuum bottoms coking or visbreaking. 
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3.2.6 Miscellaneous Stream Quality Factors 


The following input variables pertain to quality factors which 
relate the properties of refinery process unit product streams to the 
properties of their input streams. They are not required as input 
unless it is desired to replace the data base values contained within 
the program. All input variables are real numbers (non-arrays). All 
input items in this section are included between input header cards 
&UDATA1 .and &END. 

The following are API gravity differentials between the indicated 
product stream and the process unit feed. 

Quality Variable Description 

DAHKG* API of gasoline minus kerosene hydrotreater feed 

(petroleum) 

API of treated kerosene minus kerosene hydrotreater 
feed (petroleum) 

API of visbreaker light gasoline minus vacuum 
bottoms feed 

API of visbreaker heavy gasoline minus vacuum - 
bottoms feed 

API of visbreaker gas oil minus vacuum bottoms feed 

API of visbreaker pitch minus vacuum bottoms feed 

API of coker light gasoline minus vacuum bottoms feed 

API of coker heavy gasoline minus vacuum bottoms feed 

API of thermal cracker gasoline minus gas oil feed 

API of thermal cracker thermal fuel minus gas oil feed 

API of catalytic cracker gasoline minus gas oil feed 

API of catalytic cracker light cycle oil minus gas oil 
feed 

DACCS API of catalytic cracker bottoms product minus gas 

gas oil feed 


DAVBLG 

DAVBHG 

DAVBGO 

DAVBPH 

DACOLG 

DACOHG 

DATCG 

DATCTF 

DACCG 

DACCL 


* 


Light Kerosene: 400 to 525°F boiling range 



Quality Variable 

Description 

DADDD** 

API of desulfurized distillate minus distillate de- 
sulfurizer feed (petroleum) 

DADDG** 

API of gasoline product minus distillate desulfurizer 
feed (petroleum) 

DAGDG 

API of gasoline product minus gas oil desulfurizer 
feed (petroleum) 

DAGDD 

API of distillate product minus gas oil desulfurizer 
feed (petroleum) 

DAGDGO 

APIof treated gas oil product minus coker gas oil 
desulfurizer deed. 

DAGGO 

API of treated gas oil product minus straight run 
gas oil desulfurizer feed (petroleum) 

DADHH 

API of heavy hydrocrackate (250 to 400°F) minus 
distillate hydrocracker feed 

DAHYLK 

API of light hydrocrackate minus gas oil hydrocracker 
feed (in gasoline mode) 

DAHYHK 

API of heavy hydrocrackate minus gas oil hydrocracker 
feed (in gasoline mode) 

DAHYML 

API of light hydrocracker minus gas oil hydrocracker 
feed (in middle distillate mode) 

DAHYMH 

API of heavy hydrocrackate minus gas oil hydrocracker 
feed (in middle distillate mode) 

DAHYKL 

API of light kerosene (400 to 525°F) minus gas oil 
hydrocracker feed (distillate operation). 

DAHYKH 

API of heavy kerosene (525 to 650°F) minus gas oil 
hydrocracker feed (distillate operation). 




** Heavy Kerosene: 525 to 650° F boiling range. 
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Shale Oil 

Gas Oil Hydrocracking 


Coal Oil 

Gas Oil Hydrocracking 


Light 

Kerosene 

Product 


(400-525°F) 


Heavy 

Kerosene 

Product 


(525-65Q°F) 


Light 

Kerosene 

Product 


(400-525°F) 


Heavy 

Kerosene 

Product 

(525-650°F) 


Paraffin Content 



PNHLKC 

PNHHKC 

(WT5S)* 

Aromatic Content 

PNHLKS 

PNHHKS 

ARHLKC 

ARHHKC 

(WT%) 

ARHLKS 

ARHHKS 

Naphthene Content 
(WT%) 

NAHLKS 

NAHHKS 

NAHLKC 

NAHHKC 

Freezing Point 
CF)* 

FRHLKS 

FRHHKS 

FRHLKC 

FRHHKC 


API Difference 
Between Product 
and Gas Oil Feed 


DACHG1 

DACHG2 

DACHG3 


Gas Oil 

- Low Severity 

- Medium Severity 

- High Severity 


DASHG1 

DASHG2 

DASHG3 


Kerosene 

- Low severity 

- Medium severity 

- High severity 


DASHL1 DASHH1 
DASHL2 DASHH2 
DASHL3 DASHH3 


DACHL1 DACHH1 
DACHL2 DACHH2 
DACHL3 DACHH3 


Light Gasoline 

(250°F EP)* DASHLG 

Heavy Gasoline 

(250-400°F) DASHHG 


DACHLG 

DACHHG 



* Make these inputs consistent with the specified severity. 
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PETROLEUM PROCESSING 


Kerosene Hydrotreater Distillate Hydrotreater 


Fractional Change in Paraffin Content 


DPNLKP 

DPNHKP 

Fractional Change in Aromatic Content 


DARLKP 

DARHKP 

Change in Freezing Point (°F) 


DFRLKP 

DFRHKP 

API difference, gasoline product 

minus feed 


DAHKG 

DADDG 

API difference treated kerosene or 
distillate minus feed 


DAHKK 

DADDD 


Petroleum Gas 

Oil Hydrocracker 


Gasoline Mode 

Distillate Mode 


— 

Opera ti on 

Operation 

Light Hydrocrackate Product 

(IBP to 250°F)* 


DAHYLK 

DAHYML 

Heavy Hydrocrackate Product 
(250 to 400°F)* 


DAHYHK 

DAHYMH 

Light Kerosene Product 

(400 to 525°F)** 

API kerosene product minus gas oil 

feed 

. 

DAHYKL 

Paraffin content, wt% 


- 

PNHLKP 

Napthene content, wt% 


- 

NAHLKP 

Aromatics content, wt% 


- 

ARHLKP 

Freezing Point, °F 


- 

FRHLKP 

Heavy kerosene product 
(525 to 650°F)** 

API kerosene product minus gas oil 

feed 

- 

DAHYKH 

Paraffin content, wt% 


- 

PNHHKP 

Napthene content, wt% 


- 

NAHHKP 

Aromatics content, wt% 


- 

ARHHKP 

Freezing Point, °F 


- 

FRHHKP 

* API gravity of gasoline range products 

minus the gas oil feed 


** Kerosene products are not yielded in gasoline mode operation. 

3.2.7 Purchased Blending Stocks and Fixed Product Sales 

The following input variables relate to the use of segregated 

i 

refinery blending stocks. Values must be input if any of these segregated 
feed stocks are used since of course no default values are contained 
within the program. These streams are used directly for blending and 
are not "processed" in the refinery. All input in this section must be 
included between the &PDATA and SEND header cards. 
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Description 

Volume (BPD) of segregated kerosene stock (400 to 525°F). 

API of segregated kerosene stock. 

Sulfur content of segregated kerosene stock, wt percent. 
Viscosity of segregated kerosene stock, viscosity blending 
index at 210°F. 

Nitrogen content of segregated kerosene stock, wt percent. 

Hydrogen content of segregated kerosene stock, wt percent. 

Paraffin content of segregated kerosene stock, wt percent. 

Naphthene content of segregated kerosene stock, wt. percent 
Aromatics content of segregated kerosene stock, wt percent. 
Freezing point of segregated kerosene stock, °F. 

Smoke point of segregated kerosene stock, millimeters. 

Heat of combustion of segregated kerosene stock, BTU per pound. 

The input logic is identical for purchased heavy kerosene or diesel 

cut (500 to 650°F), purchased gas oil (650 to 1050°F) and purchased residuum 

(1050°F +). The corresponding array names are, respectively: 

PDSL (I) 

P60IL (I) 

PRSD (I) 


There is provision in the program for fixed product sales of light 
and heavy kerosene (diesel) and gas oil. The symbolic variable names are given 


below. 




Name 

Ml 

Units 

Description 

KERO 

Real 

BPD 

The external sale of the refinery light kerosene stream. 
If there is a kerosene hydrotreater, it is assumed 
that the hydrotreated kerosene stream is being sold. 
■Boiling range is 400 to 525°F. 

DIESEL 

Real 

BPD 

The external sale of the refinery diesel oil (heavy 
kerosene) stream. If there is a middle distillate desul 
furizer, it is assumed that the desulfurized diesel 
stream is sold. Boiling range is 525 to 650°F. 

GASOIL 

Real 

BPD 

The external sale of the refinery vacuum gas oil stream 
(untreated). Boiling range is 650 to 1050OF. 


Variable 

PKER0 (1) 
PKER0 (2) 
PKER0 (3) 
PKER0 (4) 

PKER0 (5) 
PKER0 (6) 
PKER0 (7) 
PKER0 (8) 
PKER0 (9) 
PKER0 (10) 
PKER0 (11) 
PKER0 (12) 


3.2.8 Aviation Turbine Fuel and Fuel Oil Blend Specifications 

The following input variables relate to the establishment of aviation- 
turbine fuel blend specifications. If the hydrogen content specification is 
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omitted, the program will bypass blend calculations for determining the 
maximum quantity of fuel available of a fixed hydrogen content. If the 
endpoint specification is omitted, a default value of 588°F will be 
assumed midway between the 525 and 650° F distillation temperature for 
which blend component yields and properties are contained in the program 
data bases. All of the variables listed below are real numbers. All 
information in this section must be included between the &PDATA and &END 
header cards. 

Blend Variable Description 

HDSPEC The aviation turbine fuel hydrogen 

content in weight percent. 

EPSPEC The ASTM D-86 endpoint in °F specified 

for the jet fuel blend (this blend 
is calculated in addition to the 525 
and 650°F endpoint' blends, which are 
calculated as a matter of course). 

The following variables relate to the establishment of middle 
distillate and residual fuel oil blend specifications. If they are not 
input, as will usually be the case if aviation turbine fuel blending is 
the prime concern, the program will simply calculate and report the 
volumes and properties of the remaining total middle distillate and 
residual fuel oil pools. The following variables must be input if 
additional fuel oil pool calculations are desired to meet any particular 
‘specifications. All of the variables are real numbers. 

i 

Blend Variable Description 

SPECMD The program will calculate and report 

the maximum volume of middle distillate 
range fuel oil of the specified sulfur . 
content available from the entire 

middle distillate fuel oil pool. 
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SPECRF 


The program will calculate and report 
the maximum volume of residual fuel oil 
of the specified sulfur content available 
from the entire residual fuel oil pool. 

SPM1, VM1 These input variables refer to the 

SPM2, VM2 sulfur contents (SPM1, SPM2) and BPD 

volumes (VM1 , VM2) of fixed middle 
distillate range fuel oil requirements. 
Either one or two fixed blends may be 
specified. The program will calculate 
the composition and properties of these 
blends and will blend the excess middle 
distillate components into the residual 
fuel oil pool. 

These are input variables corresponding 
to the above, but which refer to the 
residual fuel oil pool. The program 
will calculate the composition and 
properties of these residual fuel oil 
blends, along with the composition, 
volume, and properties of the remaining 
residual fuel oil pool after the specified 
fixed blends have been subtracted out. 

3.2.9 Report Options and Economic Data 

The following report options are available for specifying the type and 

quantity of output reports: 

IREP Integer IREP=1, Sunmary reports and fuel blending 

reports only. This option excludes the 
refinery economic summary report. This 
is the default option; IREP=2, the above 
reports, plus the refinery economic 
summary report; IREP=3, all reports, 
including material and sulfur, nitrogen, 
hydrogen balances for each of the refinery 
process units. 

. Economic Inputs 

These are required only if report options IREP=2 or 3 are in force. 

Each of the input items entering into the economic calculation are specified 
below. Required items are listed first, followed by those which are optional 


SPR1 , VR1 
SPR2, VR2 
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because of the inclusion of corresponding values within the program data 
base. Economics related data base values are given in Exhibit 7. All 
information in this section must be included between the &PDATA and SEND 
header cards. 

Required Inputs 
Symbol i c 


Name 

1m. 

Units 

Description 

STRF 

Real 

Number 

Fraction 

Average refinery stream day factor. 
Fraction of time that the refinery 
is not down for maintenance or 
emergency repair. 

PDBBL (1) 

Real 

dollars per 

Refinery netback price for fuel gas 


Array 

barrel fuel 
oil equivalent 
(FOE) 


PDBBL (2) 

1! 

dollars per 
liquid barrel 
@ 60°F 

Refinery netback price for LPG sales 


; 
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Symbolic 

Name 

Im. 

Units 

PDBBL (3) 

Real 

Array 

dollars per 
liquid barrel 
0 60°F. 

PDBBL (4) 

II 

II 

PDBBL (5) 

fl 

11 

PDBBL (6) 

II 

It 

PDBBL (7) 

II 

II 

PDBBL (8) 

II 

11 

PDBBL (9) 

II 

II 

PDBBL (10) 

II 

II 

DPBBL (I) 

II 

dollar per 
liquid barrel 
0 60°F 

PFPR 

real 

number 

dollars per 
barrel of 
fuel oil 
equivalent 
(FOE) 

CENTKW 

Real 

number 

cents per 

kilowatt 

hour. 

UNTMP (I) 

II 

Percent 


Description 

Refinery netback price for naptha export. 


Refinery netback price for refinery 
kerosene sale (if any). 

Refinery netback price for refinery 
diesel sale (if any). 

Refinery netback price for refinery 
gas oil sale (if any). 

Refinery netback price for the total 
gasoline pool produced. 

Refinery netback price for the specified 
aviation turbine fuel blend. 

Refinery netback price for middle 
distillate fuel oil. 

Refinery netback for residual fuel oil. 

Crude oil cost, where index I refers 
to the crude oil to be given in Exhibit 
4, List of Data Base Crude Oils. If 
crudes outside of the data base are 
specified, the codes should correspond 
to those used for the newly specified 
crude oils. 

The price assigned to fuel required to meet 
the internal refinery fuel needs. If these 
are to be diverted from potential refinery 
fuel oil or gas sales, then the refinery 
netback price for these may be used. If 
outside fuel purchases are required, then 
the purchase price may apply. 

The refinery power purchase cost. 


The percentage of investment cost to 
be taken for fixed annual maintenance costs 
for the specified process unit. The 
index I refers to the process unit 

identifier given in Exhibit 5 and is 
used to specify the unit. 
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dollars per Process unit variable operating costs, 
barrel of feed 


dollars Process unit fixed labor costs 

per day. 

Fraction The annualized refinery investment carryinc 

charge, expressed on a fraction rather 
than a percentage. This carrying charge 
usually includes the cost of capital, local 
taxes, and plant depreciation. 


Optional Inputs (Values are contained in the program data base) 


Symbolic 

Name 

1m. 

Units 

Description 

UNTCPC (I)* 

Real 

array 

dollars per 
barrel per 
day ($ per 
BPD of feed) 

Process unit capacity cost corresponding 
to the unit size, input below. 

UNTSTR (I)* 

ll 

BPD 

The unit size for which the process 
unit cost is given. 

UNTEXP (I) 

11 

dimensionless 

fraction 

The power law exponent to be used in 
adjusting the process unit cost for 
actual plant size. The equation is: 



I is investment cost 
C is BPD capacity and 
N is the power law exponent 


UNTVC (I)* 

UNTLC (I)* 

IVESTC Real 

number 


Exceptions are the alkylation unit, which is based on $/BPD of 
alkylate product, the polymerization unit which is based on $/BPD of 
polymer product and the hydrogen plant which is based on dollars per 
thousands of standard cubic feet per day (60° and 14.696 psia). Note 
also that the index I refers to the process unit identifier given in 
Exhibit 5. 
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3.3 INPUT VARIABLE DEFAULT VALUES 


Input variables are given default values within the program in 
order to minimize the data input requirements. Please note that a 
subject case of a case study sequence will use the data input established 
in the base case, rather than revert to the original default values. 

For example, if a base case has a fluid catalytic cracker, FCC = 50000.0, 
and it is desired to study the refinery situation without the cracker, 
then FCC = 0.0 must be input in the subject case, rather than depending 
on the default value of FCC = 0.0. In brief, the base case values 
override the program default values. 

3.3.1 Unit Throughputs 

All unit throughputs are initialized to 0.0, so that if a unit is 
not present within a refinery configuration, no input item is required 
for that unit. The atmospheric crude unit capacity is calculated by the 
program to accomodate the total quantity of crude oil charged to the 
refinery. Vacuum distillation capacity, however, must be specified. 

3.3.2 Unit Operations 

a) CONV = CONVH =85.0 

REFRON = RONHI = 100.0 

Fluid catalytic cracker and catalytic Naphtha Reformer severities 
are initialized at the maximum. 

b) MODE = 1 

The Gas Oil Hydrocracker is assumed to be in the gasoline production 
mode unless otherwise specified. 

c) SVSKHT, SVCKHT, SVSGHT, SVCGHT all = 2. 
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The Shale Oil and Coal Oil Kerosene Hydrotreater and Gas Oil 
Hydrocrackers are all initialized at medium severity. 

3.3.3 Internal Stream Destinations 

The default option is to prorate the volume of internal refinery 

streams having more than one downstream destination by the capacity of 

the alternative downstream process units. For example, if: 

GODES = 50000.0 Gas Oil Desulfurization Capacity 

FCC = 50000.0 Fluid Catalytic Cracker 

then: 

VGODES =0.5 Fraction virgin gas oil to desulfurization 
CGODES =0.5 Fraction coker gas oil to desulfurization 
VGOFCC =0.5 Fraction virgin gas oil to cracker 
CGOFCC =0.5 Fraction coker gas oil to cracker 

In the above, any gas oils which can not be accomodated by the 

process units because of restrictive capacity will automatically go to 

residual fuel oil blending. If the fractional disposition of a given 

internal refinery is less than 1.0, then the remainder will automatically 

be sent to final product blending, regardless of process unit capacities. 

For example, if the only destination for light cycle oil is desulfurization, 

but it is desired to send 20% directly to middle distillate fuel oil 

blending, then input LCODES =0.8 and the program automatically sends the 

remainder to fuel oil blending. 

3.3.4 Product Sales 
KERO =0.0 
DIESEL = 0.0 
GASOIL =0.0 

It is assumed that unless otherwise specified there are no external 
kerosene, diesel or gasoil sales. 
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3.3.5 Segregated Blending Stocks 


PKERO (1) = 0.0 
PDSL (1) » 0.0 
PGOIL (1) = 0.0 
PRSD (1) = 0.0 

The volumes of purchased segregated fuel oil blending stocks (which 
remain unmixed with any other refinery streams prior to blending) are 
initialized at zero. The actual volumes and stream properties must be 
supplied by the program user if segregated blending stocks are desired. 
3.3.6 Blend Specifications 

The following blend specification values are present for aviation 
turbine fuel-: 

EPSPEC = 588. 0°F 
HDSPEC =0.0 


The endpoint specification default is set midway between 525 and 
650°F, the ASTM temperature corresponding to blending component pro- 
perties contained within the program data base. The effect of the 
hydrogen specification default of 0.0 is to bypass program options that 
blend aviation turbine fuels to a minimum hydrogen content blend specification. 

The following blend specifications values are present for middle 
distillate and residual fuel oil blends: 


SPECMD =0.0 
SPM1 =0.0 
SPM2 =0.0 
SPECRF =0.0 
SPR1 = 0.0 
SPR2 =0.0 


VM1 =0.0 
VM2 =0.0 

VR1 =0.0 
VR2 =0.0 
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The effect of the above zero Initialization values is to make the 
programmed assumption that no fixed volume and sulfur content blends are 
to be produced, unless explicitly specified by the program user by using 
one or more of the above input variables. If none of the above input 
variables are specified, the program will simply calculate the volume 
and properties of the distillate and residual fuel oil pools, after jet 
fuel production. 

3.3.7 Data Base Values 

Yields, stream properties, and stream quality factors are contained 
within the program to approximate a typical petroleum refinery with a 
wide range of individual process units. In addition, crude assay data 
is contained within the program for a number of crude oils of world-wide 
importance. The full list of stored data is not presented here. The 
complete data base is contained in the programming documentation found 
in Volume III (NASA CR-1 35335). 

Any item in the process unit and crude assay data base may be 
altered by the user to represent variations in process unit yields, 
stream qualities, and crude oil assays. This affords the model user 
considerable power in matching a particular petroleum refining situation. 

Crude oil and process unit data base items must be changed at the 
beginning of the input data stream run set up and remains in force for 
all subsequent base and subject cases throughout the run sequence. 

Separate job steps are required if it is desired to study the 
effect of varying the crude assay or process unit data base values. The 
crude oil and process unit data base files are handled by the program as 
permanent files throughout a run sequence once the initial changes are 
made at the beginning of the input card stream. 
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Process unit investment costs, fixed and variable operating cost 
information, and investment cost power law exponents are contained In 
the program data base (in subroutine Block Data). These values are used 
in the economic calculation, unless altered by the user during run 
sequence. 
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4.0 MODEL FLEXIBILITY AND LIMITATIONS 


The Gordian Refinery Simulation Model was developed using the 
Control Data 6600 computer and has been converted to the IBM 360/67 
computer. Their program requires about 225 K words of storage on the 
CDC 6600 computer and about 400 K bytes of storage on the IBM 360/67 
computer. 

The detailed model limitations are covered either explicitly or 
implicitly in section 3.0 (the input description) of this manual. In 
general, any capability not encompassed by the description of program 
inputs is not available to the model. The purpose of this section is 
to summarize the principal model limitations and flexibility aspects. 

(1) The total crude data base is limited to 35 distinct crude 
types. The program is presently organized that 26 of these are reserved 
for permanent data base crudes and 9 are for crudes which the user may 
want to add during a particular run sequence. Petroleum based, coal 
based and shale oil based crude oils may be specified for refinery pro- 
cessing. 

The word crude is used somewhat loosely. Gas oils, diesel stocks 
and kerosenes may also be made a part of either the permanent or temporary 
data base. The capability reserved for additional crude types may also 
be used to represent refinery purchase streams such as natural gas and 
butanes. 

(2) An important factor in using the above refinery feed capab'ilites 
is to recognize that petroleum, coal based and shale oil based crude oils 
are processed through segregated crude units and hydrotreating facilities. 

In order to simulate the situation where segregated blending stocks are used, 
as many as four purchased blending stocks may be specified. 
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The input terminology suggests that the purchased streams are the 
kerosene, diesel, gas oil, and residuum types. Actually, since all of 
the stream properties must be supplied along with the volume rate of 
the feed, these streams may represent any segregated stream in general. 

Thus, it would be possible to specify two segregated gas oil and two 
segregated kerosene streams as purchase for final fuel oil blending stocks. 

Any streams of the above type are not blended with any other feeds 
prior to being blended into their respective fuel oil pools. 

(3) The refinery process units which are included in the model's 
process unit data base are restricted to those specified in Exhibit 5. 
Programming modifications are necessary to expand the allowable 
petroleum refinery configuration beyond its present range. 

The model is flexible in allowing for diversity in the specification 
of the refinery configuration. For example, units using the same feed 
stream type such as a fluid cracker and a thermal cracker (using straight 
run virgin gas oil) may be included within the same refinery configuration. 
When dual units of this type are included, it is necessary to define the 
split in the feed stream between units by making use of the appropriate 
input variables, unless the default of prorating feeds by available capacity 
is desired. 

(4) Aviation Turbine fuel (jet fuel blends) in the kerosene range 
are reported by the model. Considerable flexibility has been included in 
the program in order to extract the maximum useful information from a 
single run. The user may specify a jet fuel endpoint specification and 
minimum hydrogen content specification. The model will then calculate 
blends as 525?F, 650?F and the specified endpoint. All blends are reported 
on an "as made" basis and at the specified hydrogen content. 
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The aviation turbine fuel of specified endpoint and specified hydrogen 
content is calculated as the basis for the overall refinery balance, with 
the remaining distillate components going to distillate fuel (heating oil) 
blending. 

(5) Up to two fixed distillate blends and two fixed residual fuel 
oil blends may be specified. Both sulfur content and blend volume must 
be specified for each fixed blend. 

Each distillate blend specified is subtracted from the total distillate 
fuel oil pool and the volume, composition, and properties of the remaining 
pool are calculated by the model. Any excess middle distillate components 
are then blended into the residual fuel oil pool. The required residual 
fuel oil blends are then made and the properties of the remaining residual 
pool are reported. 
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5.0 SAMPLE PROBLEM SET-UPS 


The copies of computer card input shown in Exhibit 8 illustrate 
the feature of the Gordian Refinery Simulation Model of handling a base 
case and subsequent dependent cases. The first group of images is the 
sample case presented in Exhibit 1. Two dependent input setups are shown 
for changes in jet fuel, distillate, and residual fuel outputs. Each 
input group represents a change only in the variables shown; all others 
remain at base case values. The third case has an input variable to 
indicate that it is the last of the dependent series and the subsequent 
run will be a new base case. 
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6.0 DIAGNOSTICS 


Input editing for system errors is accomplished by the FORTRAN 
Namelist feature. Such errors as mispunched header cards, mispunching 
or misspelling of input variable names, omission of commas and errors in 
numeric field punching are flagged by Namelist and generally result in 
problem termination. An additional diagnostic has been included in the 
program to signal a premature end-of-file in the input stream, i.e. 
before an SEND trailer card has been reached. 

The program has been written to provide automatically for such 
contingencies as undersized process units, and deficiencies in hydrogen 
or fuel availability by blending off excess feed streams and by indicating 
the need for purchase. The greatest propensity of the program to fail 
in the calculation is probably in the jet fuel, distillate and residual 
fuel oil blending calculations. The problem arises in trying to make 
specified fuels outside of the ability of the program to do so — for 
example, if a hydrogen specification content of 15.0 weight percent were 
specified for a 650°F endpoint jet fuel blend and the highest blending 
component hydrogen content were 14.8 weight percent. 

The subroutine involved in making the blending calculation is BMAX. 
Either or both of two messages are issued by BMAX if a blending problem 
develops: 

"A FUEL OIL BLEND CALCULATION HAS FAILED — CHECK INPUT" 

"A ZERO DIVIDE HAS OCCURED IN SUBROUTINE BMAX — THE RESULTS ARE 

THEREFORE SUSPECT — CHECK PROBLEM INPUT" 

The input checks under this circumstance would usually consist of 
redefining a feasible problem, i.e., a refinery configuration capable of 
making the specified blends. 
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7.0 CONCLUDING REMARKS 


This report has described the Gordian Refinery Simulation Model , 
a Fortran computing program for predicting the flow streams and material, 
energy, and economic balances of a typical petroleum refinery, with 
particular emphasis on production of aviation turbine fuel of varying 
end point and hydrogen content specifications. Program logic, optional 
and required input, default values, diagnostics, and sample output are 
discussed. A brief example is included, illustrating the case feature 
for parametric or optimization study variations of a base case. 

It is felt that the material presented in this report will be 
sufficient for the needs of many readers. The simplicity and flexibility 
of the program can be readily assessed, and the limitations should be 
apparent, particularly with respect to uncertain reliability when dealing with 
the properties of shale and coal syncrudes. The sample input and output 
printouts presented here can serve as a guide for interpretation of future 
reported listings and data. Additional details, however, will be found 
in Volume II of this report series (NASA CR-135334) covering the mathe- 
matics and data bases and Volume III covering the program listing, both 
documents to be available from the NASA Project Manager. Computer docu- 
mentation and tapes can be purchased through the Computer Software 
and Management Information Office (COSMIC), 112 Barrow Hall, University 
of Georgia, Athens, GA, 30602, under the number LEW-13047. 
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EXHIBIT 1 


COMPUTER PRINTOUT FOR A SAMPLE CASE 
On the following pages, a printout of a sample case using the 
Gordian Refinery Simulation Model is reproduced. On the first page, the 
optional input data are illustrated. A crude assay data change input is 
shown between &CDATA and &END, an alkylation unit feed nitrogen content 
change between &UDATA and SEND, and an; alkylation unit product octane 
number change between &UDATA1 and &END. (The printing unit lacks a 
& symbol and substitutes a $ on the printout). On the second page, 
required data between &PDATA and SEND are shown. The selected case 
represents a refinery model with crude feed of 33,000 BPD each of Agha 
Jari petroleum, Paraho shale, and Kentucky synthoil coal oil. The 
input variables prescribe the sizes of the process units, catalytic 
cracker conversion, mid-distillate operation, and medium hydrotreater 
and hydrocracker severity. These input values are summarized on a page 
titled, "SUMMARY OF REFINERY UNIT OPERATIONS" near the end of the print- 
out pages. The input IREP=3 specifies a complete printout, including 
all material balances and economic calculations. The production of jet 
fuel pools with a specified 600® F end point and 12.8% hydrogen is indicated. 
This sample case is discussed briefly in Section 2.2 of this report. 
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EXHIBIT 2 - Block Flow Diagram of Program Logic 
MAIN PROGRAM 








UNTENR 
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ECON 


Calculate the 
oil refinery 
investment and 
operating costs 


Generate 
summary reports 


I 


Return to the 
beginning of the 
main program for 
successive cases in 
the run sequence 


1 


Calculate the overall 
refinery material 
balance, energy usage 
and gasoline blending. 

Generate summary 
material balance, 
energy efficiency, 
refinery operation 
and gasoline blending 
summary reports 


T 


BLENDF 

(cont.) 


Call auxiliary* 
subroutines 


BLDREP 



Generate 
blend reports 


*BLDARY - Sets up program arrays 
containing blending stream volumes 
and properties and calculates the 
properties of jet fuel blending 
components of varying endpoints 

SSORT & SSORTI - Sorts blending 
components by sulfur and 
hydrogen content ranking 

BMAX - Calculate maximum 
blend quantities of specified 
hydrogen and sulfur content 

BPROP - Calculate all' jet fuel, 
distillate and fuel oil blends 


EXIT at the end of the 
last case 
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EXHIBIT 3 


REFINERY SIMULATION BLOCK FLOW DIAGRAM 



FF.-rrvrrrr-r |j PVrrrrrrrrr 

























EXHIBIT 4 


LIST OF DATA BASE CRUDE OILS 


PROGRAM INPUT 
IDENTIFICATION 

CODE CRUDE OIL 


1 

Tigre-Venezuelan - 

1 .57%S, 

24.7 API 

2 

Lot 17-Venzuelan - 

0.98%S, 

36.1 API 

3 

Bachequero-Venezuel an 

2.4%S, 

16.8 API 

4 

Nigerian Light - 

0.1 4%S 

34.7 API 

5 

Amal -Libyan - 

0.10%S, 

35.8 API 

6 

Arzew-Algerian - 

0.10%S, 

44.1 API 

7 

Bakr-Egyptian - 

4.4%S, 

19.6 API 

8 

Light Arabian - 

1 .65%S, 

34.2 API 

9 

Agha Jari- Irani an - 

1 . 34%S , 

34.3 API 

10 

Kuwait - 

2.53%S, 

31.4 API 

11 

Paraho Shale Oil - 

0.7US, 

19.3 API 

12 

Tosco Shale Oil - 

0.67%S, 

21.0 API 

13 

Garrett Shale Oil 

0.64%S, 

25.0 API 

14 

Kentucky Coal Synthoil - 

0.22%S, 

5.9 API 

15 

Prudoe Bay Alaska - 

1 . 04%S , 

26.8 API 

16 

Ekofisk-North Sea- 




Norwegian - 

0.18%S, 

35.6 API 

17 

West Texas Sour - 

1 . 9%S , 

34.0 API 

18 

South Lousiiania 




Ostrica - 

0.31%S, 

32.3 API 

19 

Louisiana Delta - 

0.30%S, 

30.6 API 

20 

East Texas 

0.3%S, 

38.0 API 

21 

Utah Aneth - 

0.12%S, 

40.9 API 

22 

Wyoming Sour - 

2.4%S, 

24.9 API 

23 

Oklahoma Golden Trend - 

0.20%S, 

39.9 API 

24* 

California Elk Hills - 

0.61%S, 

26.6 API 

25 

California Wilnington - 

1 .43%S, 

21.7 API 

26 

Canadian Pembina - 

0.83%S, 

32.7 API 

27-35 

Assay Data, as input 




for additional crude 
oils 


i 


* Incomplete assay data - do not use without inputing assay data. 
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EXHIBIT 5 

INTERNAL REFINERY STREAM AND PROCESS IDENTIFIERS 


Program 

Identification 

Number • Description 


1 

2 

3 

4 

5 

6 

7 

8 
9 

10 

11 

12 

13 

14 

15 

16 

17 

18 

19 

20 
21 
22 

23 

24 

25 

26 

27 

28 

29 

30 

31 

32 

33 

34 

35 

36 

37 


Refinery Steam Consumption 
Refinery Fuel Consumption 
Refinery Power Usage 
Hydrogen (BFOE) 

Ethane and lighter (BFOE) 

Propane 

Propylene 

Isobutane 

n-Butane 

Butylene 

Pentanes (from crude oil) 

Reformer Feed Preparation Overhead (C5/C6) 
Distillate Hydrocracker Pentanes (IC5/NC5) 
Heavy Hydrocrackate (from both distillate 
and gas oil ) 

Light Hydrocrackate (from gas oil hydro- 
cracker) 

Reformer Feed Preparation Bottoms 
Distillate product from gas oil hydro- 
cracker 

Light Reformate 
Heavy Reformate 
Total Reformate 
Light Alkalate 
Heavy Alkylate 
Polymer Gasoline 
Light Straight Run Gasoline 
Heavy Straight Run Naphtha 
Petroleum Light Kerosene 
Petroleum Heavy Kerosene 
Petroleum Vacuum Gas Oil 
Vacuum Bottoms Reduced Crude 
Reduced Crude 

Light Coker/Visbreaker Gasoline 

Heavy Coker/Visbreaker Gasoline ' 

Light Coker/Visbreaker Gas Oil 

Heavy Coker/Visbreaker Gas Oil 

Total Coker/Visbreaker Gas Oil 

Coke produced from Coker 

Full range cracked gasoline (FCC/TCC) 


73 



38: 

39 

40 

41 

42 

43 i 

44, 

45' 

46 


47 


48 

49 

50 

51 ; 

52 

53 

54 

55 

56 

57 

58 

59 

60 
61 
62 

63 

64 

65 

66 


Light cycle oil from gas oil cat cracker 
Tower bottoms from gas oil cat cracker 
Desulfurized light cycle oil 
Hydrotreated straight run diesel 
Desulfurized light coker gas oil 
Distillate product from Gas Oil 
desulfurizer 

Gas oil product from Gas Oil 
desul furizer 

Gasoline byproduct from all 
desul furization 

Stream 43, but with Coker/Visbreaker 
Feed Source 

Stream 44, but with Coker/Visbreaker 
Feed Source 

Stream 17, but with Coker/Visbreaker 
Feed Source 

Composite Petroleum Crude Oil 
Hydrotreated Petroleum Light Kerosene 
Visbreaker Pitch . 

Thermal fuel produced from the Gas Oil 

Thermal cracker (all 400+ material) 

Composite Shale Oil 

Composite Coal Syn Crude Oil 

Shale Oil Light Kerosene 

Shale Oil Heavy Kerosene 

Shale Oil Gas Oil 

Hydrotreated Shale Oil Light Kerosene 
Hydrotreated Shale Oil Heavy Kerosene 
Hydrotreated Shale Oil Gas Oil 
Coal Oil Light Kerosene 
Coal Oil Heavy Kerosene 
Coal Oil Gas Oil 

Hydrotreated Coal Oil Light Kerosene 
Hydrotreated Coal Oil Heavy Kerosene 
Hydrotreated Coal Oil Gas Oil 


t 
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EXHIBIT 5 (cont'd) 


INTERNAL PROCESS UNIT IDENTIFIERS 


Identi fi cation 
Number • 


Description 


1 

2 

3 : 

4 

5 

6 

7 

8 
9 

10 

11 

12 

13 

14 

15 

16 

17 

18 

19 

20 
21 
22 


Total crude unit capacity 
Visbreaker 

Catalytic Cracker (petroleum) 
Thermal Cracker (petroleum) 
Alkylation Unit 
Polymerization Unit 
Distillate Desulfurizer (petroleum) 
Gas Oil Desulfurizer (petroleum) 
Distillate Hydrocracker (petroleum) 
Gas Oil Hydrocracker (petroleum) 
Catalytic Reformer 
Hydrogen Plant 
Butane Isomerization Unit 
Fluid Coker 

Kerosene hydrotreater (petroleum) 
Crude Unit (petroleum feed) 

Crude Unit (shale oil feed) 

Crude Unit (coal oil feed) 

Shale Oil Hydrocracker (gas oil) 
Coal Oil Hydrocracker (gas oil) 
Shale Oil Hydrotreater (kerosene) 
Coal Oil Hydrotreater (kerosene) 


EXHIBIT 6 

Linear Viscosity Blending Chart 


NOTE: Use this chart at 210°F fcr 

use with the Gordian Fuels 
Blending Model. Conversions 
may be made from either Saybolt 
Seconds Universal (SSU) or 
Saybolt Seconds turol (SSF) 
viscosities. 


3 
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EXHIBIT 6 (cont'd) 






TO 

25 

30 

35 

UO 

50 

60 

80 

100 

150 

200 

300 

1»50 

6oo 


1,000 


2,000 


3,500 

5,000 

10,000 


20,000 

. 35,OCO 
50,000 

100,000 


300.000 

500.000 
1,000,000 


3.000. 000 

5 . 000 . 000 
’10,000,000 


100,000, coo 


1,00^.000,000 


10 , 000 . 000.000 




250 

300 

l*C0 

500 


750 

1.000 

1,500 

2,000 

3.000 

fe.ooo 

6.000 

8.000 
10,000 


20,000 


35.000 

50.000 


100.000 

200.000 

350.000 

500.000 

1 , 000,000 

2.5CC.OOO 

5 , 000.000 

10 . 000.000 

30.000. ccc 

50. 000. COO 
lOO.OCO.CCO 
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EXHIBIT 7 


ECONOMIC DATA BASE 


1. CAPITAL INVESTMENT COSTS 

The unit capital investment costs (Jan.l, 1977 basis) stored in 


the economic subroutine are shown below: 

< 



Unit Capital 

Reference 

Refinery Unit 

Investment Cost 

Unit Size 


($/bpsd) 

(bpsd) 


Source 


Atmospheric Distillation 
Vacuum Distillation 
Catalytic Gas Oil Cracker 
Thermal Gas Oil Cracker 
Kerosene Hydtrotreater 
Gas Oil Oesulfurizer 
Distillate Oesulfurizer 
Fluid Coking 

Vacuum Bottoms Visbreaker 
Distillate Hydrocracker 
Gas Oil Hydrocracker 
Catalytic Naphtha Reformer 
Alkylation 
Polymerization 
Isomerization 
Hydrogen Manufacture 


190 

198 

699 

185 

122 

614 

267 

770 

185 

747 

904 

543 

1193 

800 


222 

170/mscfd 


20,000 

20,000 

40.000 

20.000 

30.000 

50.000 

30.000 
20,000 
. 20,000 

15.000 

20.000 
20,000 
10,000 

500 

6,650 

50 mmscfd . 


OGJ, Feb. 25, 1974, p.71 
OGJ, Mar. 4, 1974, p.100 
OGJ, Apr. 15, 1974, p.66 
OGJ, Apr. 8, 1974, p. 74 
HP, September, 1976 
OGJ, Mar.l , 1976 p. 120 
HP, September, 1976 
OGJ, May 24, 1976 
OGJ, Apr. 8, 1974 
HP, September, 1974 
OGJ, Mar. 25, 1974 p. 120 
OGJ, Apr. 22, 1974 p. 130 
OGJ, Apr. 8, 1974 p. 74 
"Gordian Estimate" 

HP, September, 1974 
OGJ, Mar. 25, 1974 p. 120 


2. ECONOMICS OF SCALE 


The following power law formula is used in the model to calculate 
economics of scale. 



Where I, I 2 are investment costs C, Co are the corresponding 
capacities in BPD and N is the power law exponent. 

Standard sizes and costs used are given in part (1) above, the 
power law exponents contained in the model appear on the next page. 
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The same exponents are assumed for processing shale and coal oil kerosene and gas oil cuts. 



3. VARIABLE OPERATING COSTS 


The catalyst and chemicals cost figures are in terms of costs per 
barrel of crude throughput. In order to update the chemicals and catalysts 
costs, the Nelson-Chemicals Index ( Oil and Gas Journal) was used. All 
chemicals and catalyst costs used in the model are on a Jan. 1, 1977 
basis. Water cost were based on an average reported consumption (gal/ 
.barrel), at a unit cost assumed to be $0.05/1000 gals. 

The total variable operating costs used in the model (chemicals, 
catalyst, and water) are shown below: 


Unit 

Atmospheric Distillation 
Vacuum Distillation 
Catalytic Gas Oil Cracker 
Thermal Gas Oil Cracker 
Kerosene Hydrotreater 
Gas Oil Desulfurizer 
Distillate Desulfurizer 
Fluid Coking 

Vacuum Bottoms Visbreaker 
Distillate Hydrocracker 
Gas Oil Hydrocracker 
Catalytic Naphtha Reformer 
Alkylation 
Polymerization 
Isomerization 
Hydrogen Manufacture 


Variable Cost ($/Bbl) 

0.015 

0.019 

0.074 

0.071 

0.086 

0.115 

0.115 

0.024 

0.022 

0.423 

0.423 

0.183 

0.990 

0.274 

0.043 

0.097 ($/mscf )' 


tHe.Tfueli power and steam consumptions are calculated, for each j unit 
based on the data base values contained in subroutine UBASE. Fuel and 
steam is produced from crude oil and intermediate streams and value is a 
function of crude oil prices. Electric power usage for the total refinery 
is calculated by the model and reported in thousands of kilowatt hours 
per day. The cost of electricity is input in cost per kilowatt hours 
by the user. 
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4. FIXED OPERATING COSTS 


The fixed operating costs include the labor and maintenance costs. 

Labor costs are assumed to be independent of throughput, units are 
assumed to be operated on a 24 hour per day basis. Labor use (men per 
shift) figures were obtained from A Guide to Refinery Operating Costs, 
by W.L. Nelson, Petroleum Publishing Company, 1976. Labor rates used 
were: operators $10.00/hr, supervisors - $12. 50/hr. In each case the 

labor rates are inclusive of benefits. 

Maintenance costs are given as a per cent of the corrected onsite 
unit capital investment cost. The computer model utilizes these maintenance 
figures to calculate the maintenance cost in dollars per day. The 
maintenance and labor cost numbers used by the computer are listed 
below: 


Unit 

Labor Cost 

Maintenance 


($/day) 

(% Capital Investment 
Per Annum 

Atmospheric Distillation 

806.10 

4.0 

Vacuum Distillation 

1002.00 

3.5 

Catalytic Gas Oil Cracker 

1545.00 

4.0 

Thermal Gas Oil Cracker 

1245.00 

5.5 

Kerosene Hydrotreater 

720.00- 

4.0 

Gas Oil Desulfurizer 

720.00 

4.0 

Distillate Desulfurizer 

720.00 

4.0 

Fluid Coking 

1485.00 

5.5 

Vacuum Bottoms Visbreaker 

1059.00 

4.3 

Distillate Hydrocracker 

1065.00 

4.0 

Gas Oil Hydrocracker 

1065.00 

4.0 • 

Catalytic Naphtha Reformer 

1012.50 

4.0 

Alkylation 

1050.00 

5.5 

Polymerization 

505.00 

4.3 

Isomerization 

1440.00 

3.5 

Hydrogen Manufacture 

585.00 

3.5 
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REFINERY OFFSITES 


Other refinery investment costs include supporting utilities, 
tankage, off limits battery piping, roads and accessways, sewers and 
drains and dock facilities. These are calculated as a function of 
Nelson refinery using the attached correlation presented by Nelson. 

The basic refinery process unit investment is multiplied by the 
percentage given to obtain the offsites portion of the total investment. 
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EXHIBIT 8 


INPUT PRINTOUT FOR SAMPLE BASE CASE AND DEPENDENT CASES 
The following page reproduces the computer printout for input 
representing the feature of base case use with subsequent dependent cases 
for parametric studies. The first group of images represent the case 
shown as Exhibit 1; the input ICASE = 1 indicates that this is a base 
case. The second group represents a subsequent calculation, in which 
the only change is a requirement of jet fuel production at 575°F end 
point and 13.0% hydrogen (compared to 600° F and 12.8% in the base case.) 

The input ICASE=2 indicates that this case is dependent on the preceding 
ICASE=1 calculation and only the indicated inputs change. The third 
group represents a second dependent calculation in which 500 BPD each 
of 0.05 and .10% sulfur middle distillate and 0.3 and .35% sulfur residual 
fuel oils are to be produced, and maximum volumes of 0.15% sulfur distillate 
and 0.6% sulfur residual are to be calculated and reported. The input 
ICASE=3 indicates that this is a dependent study of the preceding ICASE=1 , 
but it terminates the series of dependent cases. 


I 
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FUELS eiEKOIhG KOOEL-GOROIAN ASSOCIATES 
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